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Abstract 

In this paper we study scattering of two-dimensional massless Dirac fermions 
by a potential that depends on a single Cartesian variable. Depending on the 
energy of the incoming particle and its angle of incidence, there are three dif- 
ferent regimes of scattering. To find the reflection and transmission coefficients 
in these regimes, we apply the Wentzel-Kramers-Brillouin (WKB), also called 
semiclassical, approximation. We use the method of comparison equations to 
extend our prediction to nearly normal incidence, where the conventional WKB 
method should be modifled due to the degeneracy of turning points. We com- 
pare our results to numerical calculations and find good agreement. 
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In this paper we present a systematic theory of potential scattering for mass- 
less Dirac fermions. Being the effective charge carriers in graphene [1] [2 |3l |4] , 
and topological insulators O El [7], these particles attracted a keen interest. 
The discovery of massless Dirac fermions in condensed matter systems stimu- 
lated the fabrication of 'artificial graphene', a material with a hexagonal lattice, 
where quantum dots [5], or molecules [5], play the role of carbon atoms. The 
electron excitations in these materials give rise to massless Dirac fermions. The 
main feature of massless Dirac fermions is chirality (as it is called for graphene) 
or helicity (for topological insulators), i.e. an additional degree of freedom that 
relates to two kinds of particles (electrons and holes) simultaneously present in 
the system. Chirality makes the behavior of massless Dirac fermions dramat- 
ically different from that of Schrodinger particles. One of the most prominent 
examples is Klein tunneling [2l[l0l[ni[12l[l3l[T4l[T5l[l6]. Due to this effect, a 
massless Dirac fermion normally incident on an electrostatic potential will be 
transmitted with unit probability. 

In this paper we consider scattering of massless Dirac fermions by quasi- 
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Figure 1: Illustration of angular scattering by a quasi onc-dimcnsional potential barrier (gray 
area). The incoming and reflected waves are shown by arrows on the left of the potential and 
the transmitted wave is shown on the right. 



one dimensional potential barriers (the corresponding potentials depend on a 
single Cartesian variable, see figure [T]). Such barriers occur for instance in 
graphene heterostructures that were fabricated in [151 IE] ■ They can also be 
intrinsic, as in the case of puddles in graphene [171 [10]. We always assume that 
the potential profile is smooth enough, so that the Wentzel-Kramers-Brillouin 
(WKB) or semiclassical approximation [T^ [H HOI 1211 1221 123] can be used. 
The latter allows us to obtain generic formulas valid for arbitrary potentials. 
Another method that can be used to study generic potentials numerically was 
suggested in [24] . 

We distinguish three different regimes of scattering and show that the mass- 
less Dirac equation is equivalent to a pair of effective Schrodinger equations 
with complex potentials. We then solve the scattering problem for each of 
these three regimes with the help of the WKB approximation. The specific 
formulation we use is the one pioneered by Zwaan |25j . and further developed 
in [211 123 IISI im 1201 123] ■ Since we do not expect that all readers are familiar 
with this technique, it is summarized in [Appendix A| 

We start from preliminary considerations based on classical mechanics, sec- 
tion [ij In section [2] we introduce the semiclassical scattering states. Then, in 
section [3j we formulate a set of simple rules that form the basis of the WKB 
method and are sufficient to solve the scattering problem for angular scatter- 
ing. Here and further on we use the term "angular scattering" for incidence 
far from both normal and tangential (see figure [I]) . In section |4j we consider 
tunneling through a barrier supporting hole states, or an n-p-n junction. Due 
to the classically allowed hole states within the barrier, one finds Fabry-Perot 
oscillations in the transmission coefficient [TOl [El [Mj , which were used to exper- 
imentally verify Klein tunneling |15| . We show that the WKB approximation 
does not accurately describe near-normal incidence on this barrier, since the 
classical turning points are nearly degenerate in this case. To circumvent this 
obstacle and to obtain a solution that is uniformly valid in the entire range 
of incidence angles, we use the technique of comparison equations, developed 
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in [551 [Ml 1301 1311 132] and summarized in Appendix B The results presented in 
this section were already published by the authors in Q^, but no proofs were 
given there. Here we present a complete and systematic treatment, including 
the detailed derivation. 

In section [5| we apply the WKB approximation to the case of above-barrier 
scattering, and in section [6] we consider tunneling through the barrier without 
hole states. In section [t] we consider an exact solution, first constructed in [33], 
which is used to analyze the case of a single monotonous n-n junction. Surpris- 
ingly, when applying the WKB approximation to this case, we find that along 
a certain path in the complex plane this situation can be reduced to the case of 
Klein tunneling. To stress the interconnection, this case is referred to as "vir- 
tual Klein tunneling" . Finally, in section [8| we compare our predictions with 
numerical calculations. 

Our main results are presented in the form of easy-to-use analytic expressions 
for reflection and transmission coefficients. 



1. Preliminary considerations: three regimes of scattering 

The wave function '5 of a massless Dirac fermion obeys the effective Dirac 
equation 

[v(T-p + u{x/l, y/l)] ^{x, y) = E^{x, y), (1) 

where v is the Fermi velocity, cr = {ax, (Jy) is the two-dimensional vector of Pauli 
matrices, p = —ihSI is the momentum operator and I is the characteristic scale 
of change of the potential. In this paper we consider a potential u that depends 
on X only. Then the separation of variables gives "^{x^y) = (x) exp{ipyy / h) , 
and we obtain 







iPy 



Px - iPy 





u(x/l) 



(2) 



Denoting the characteristic value oi\u — E\ as upo and introducing the dimen- 
sionless variables x — x/l, Px = —ihd/dx, py — Py/pa, h = h/pol, u — u/vpo 
and E — E/vpo, we can write ([2| in the form 





Px + iPy 



Px - Wy 





- u{x) 



(3) 



or, equivalently, 

[a ■ p + u{x))-^ ^ E-^ . (4) 

Here and further on we omit the tildes. 

Let us consider the classically different scattering regimes comprised in equa- 
tion Q . It is well known [311 [H] that the classical Hamiltonian functions cor- 
responding to the matrix quantum Hamiltonian are given by the eigenvalues of 
this matrix, where momentum operators are replaced by c-numbers and correc- 
tions of the order h are neglected. Applying this prescription to Q we obtain 
two Hamiltonian functions 



iPx,x) = ±\p\ +u{x), 



(5) 
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Figure 2: Effective potentials and phase portraits in the combined representation for the cases 
E < 0: a, d; E = 0: b, e and E > 0: c, f. In the figure u(x) = —x^. 



where Lq and Lq give rise to the electron and hole dynamical systems respect- 
ively. These Hamiltonian functions coincide at the point xq if for a certain 
energy u{xo) = E. This implies Px = Py = 0, i.e. the electron and hole systems 
merge for the case of normal incidence for sub-barrier scattering. The intersec- 
tion of classical Hamiltonian functions is the origin of the Klein paradox [TOl [14] . 
It implies that the electron and hole systems cannot be treated separately for 
near-normal incidence. This forces us to change the representation to the one 
where electrons and holes are treated together. Such a representation can be 
easily found from ([5]). Indeed, for a given energy E we have 



or tIpI — 1^(2;), where we introduced the short-hand notation v{x) — u{x) — E. 
Squaring the last equality, we find 



where C{px,x) can be treated as the new Hamiltonian function and the para- 
meter e = —py plays the role of energy, so that the level lines of C{px,x) corres- 
ponding to e coincide with level lines of Lq {px , x) corresponding to the energy 
E. In the representation given by £, electrons and holes are treated together. 

The phase portraits of the Hamiltonian systems that originate from C are 
cuts of the original four-dimensional phase space \^px^Py^ x, y} by the hyperplane 
E = const. Every individual trajectory in this cut is defined by a certain 
value —Py- In figure [2] one sees the effective potentials —v'^{x) for different 
values of E and the corresponding phase portraits. These pictures describe all 



E = ±\p\+u{x), 



(6) 



^{Px,x) =pI -v'^{x) = ~pI 



(7) 



qualitatively different regimes for a Dirac particle scattered by a single hump 
potential. 

When the energy E does not exceed uq, the maximal value of u{x) (fig- 
ure [2]a,d), there exist either four (for small \py\) or two real turning points (for 
larger \py\). In the opposite case (figure [2 c,f), when E is larger than ug, real 
turning points are absent for small \py \ and there appear two of them for larger 
values of \py\. Thus we differentiate three different scattering regimes: 

1. E < uq, \py\ < uq — E: Klein tunneling regime, or tunneling through a 
barrier supporting hole states 

2. E > Uq, \py \ < E ~ Uq: abovc-barricr scattering 

3. E < Uq and \py\ > uq — i?, or > uq, \py\ > E — uq: conventional 
tunneling regime, tunneling through a barrier without hole states. 

For each of these scattering regimes we will construct a separate description. 

Representation ^ allowed us to combine electrons and holes within a single 
dynamical system. This transformation has a straightforward quantum ana- 
logue. Indeed, let us rewrite Q as (cr ■p + v{x))"^ — 0. We can now act on this 
equation from the left with the operator {a ■ p v{x)) to obtain |14| . 

(cr ■ p ~ v{x)) {(T ■ p + v{x))'ii ^ {pi +pI- v'^{x) - iha^v'ix))-^ = 0. (8) 

Since the last equation contains only a single Pauli matrix, it can be diagonalized 
by writing 

* = ( 1 \ )v2, (9) 



and one obtains 



+ v^x) -pI± zhv'ix)^ Tiia = 0. 



(10) 



The functions 771 and r/2 are not independent, they are related by: 



r]2,i ^ \h-^±iv{x)]rii^2, (H) 



Py 



as can be found from equation Q. The real part of equation (10) corresponds 
to Q, and the imaginary part gives a quantum correction to the classical trans- 
formation. 



2. Semiclassical scattering states 



Before we can solve the scattering problem for the different regimes outlined 
in the previous section, we first have to define the asymptotic scattering states. 
From now on we assume that v{z) is an analytic function in the complex plane. 



Therefore we can consider equations (10 1 and (11) for 771 in the complex plane: 



dz2 



v'^{z) -pl+ ihv'{z)^ ?7i(z) 



0, 



(12) 
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and 



Py 



h— + iv{z) ] r]i{z). 



(13) 



The semiclassical solution for equation ( |12[ ) has the form 

77i(z) = A(z,/i)e^^W/\ (14) 

where A{z, h) = Aq{z) + hAi{z) + . . . is a power series in h. Substituting this 
function into equation (12 1, we find 



h—+is'{z)^ + v\z) - + ihv' {z) 



A{z,h) = 0. 



(15) 



Equating the terms on the left hand side to zero for all powers of h, we obtain 
equations that determine s{z) and A(z, h). The terms of order h^^ give 



{s'{z)f = vHz)-pl, 
whilst collecting the terms of order h^, we find 

2s'iz)A'o{z) + s"{z)Ao{z) + v'iz)AQiz) = 0. 
Multiplying the latter equation by ^0(2;), we obtain 
d 



dz 



[s'{z)Al{z)]+v'{z)Al{z)^0. 



Assuming that s'(z) does not vanish we get 

B 



s'iz) 



exp 



.0 ^^'(0 



(16) 
(17) 

(18) 
(19) 



where _B is a constant, zq is an (up to now) arbitrary point and the integration 
should be performed along a suitable path in the complex plane. Equation ( 16 1 
has two solutions, namely ±s{zq,z), where 



s{zo,z):^±j Px{C}dC, Px{z)^ {v^{z)~pI) 



(20) 



Note that the square root is not a single-valued function in the complex plane, 
which means that we have to insert branch cuts emanating from every point 
where its argument vanishes. To distinguish the square root as an analytic 
function defined as discussed above from the positive square root of a positive 
number, we will denote the former by z^/'^^ and the latter by ^/x. Combining 
equations ( 19 ) and ( 20 ) and choosing the constants B± in an appropriate way, 
we obtain two asymptotic solutions, 



1 



exp =F 



1 



zo PxiQ 



exp ( ±-s{zq,z) 



(21) 
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The integral in the exponent can also be computed explicitly, 



dc 



v{z) 



dv 



= In 



which gives rise to the representation 
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Tl/2 



P^^iz) 



,±is{zo,z)/h 



v{z) + (v^z) ~ piy/^' 

\Py\ 



9{z) 



\Pv 



-const, (22) 



(23) 



It is important to note that g{z) does not vanish at any point z if does not 
vanish and \z\ < oo. 

At this point, let us come back to the real axis and introduce refection and 
transmission coefficients. First we establish the current conservation condition. 
It is convenient to start from equation Q, and to multiply it from the left by 
This gives 

*^(cr •p + u(a;))*(a;) = E'^\x)^[x). (24) 
Subtracting from equation (24) its complex conjugate, we obtain 

- ih^'^{x)ax^'{x) - ih[^\x)]'ax^{x) = -ihl^"^ {x)a,x^{x)]' = 0. (25) 

Hence the conserved current is given by jx{x) = {x)ax'^{x), which can also 
be written as Re [^l{x)il>2{x)] = const. Using relation (|9| between 5* and 771^2, 
we find 

hi(a:)p-h2(a:)|2 = const. (26) 



It is useful to understand the conservation equation (|26| from the point of view 

(27) 



of the effective equation (10), and equation (11). For the equation 
d^ 



/i^^ + v'^{x)-pI + ihv'{x)j 771 (x) = 0, 
the conserved quantity is given by the Wronskian, 

W 



■qi{x) Tjix) 
i][{x) r]'{x) 



(28) 



where 771(2;) and ri{x) are linear independent solutions of (27). Since, after 
complex conjugation, equation (10) for 772 coincides with equation (27), one can 
choose ri{x) = 772(0;), where the star denotes the complex conjugation. From 
([28l) we then find 



Vi{x){V2y{x) - V2{xWi{x) = const. 



(29) 



Using equation (11) to eliminate the derivatives, we arrive at (26). 



Now we introduce the scattering solutions in a given classically allowed re- 



gion. 



mix) 



ai 



\/pJ^VG{x) 



^iS(xo,x)/h ^ 



\fpJx) 



'iS{xQ .x) I h 



(30) 
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Note that we now use the symbol y^, i.e. we assume that we are on the real 
axis. We denoted 



G(x) 



\v{x)\ +Pa:{x)\'' 



\Py\ 



, V = sgn[i;(a;o)] 



and the momentum Px{x) and the action S{xq, x) are defined through 



(31) 



(32) 



In the electron region one has E > hence v(x) < 0. Applying equation (111 
to the expansion ( 30 ) , we can find r]2 {x) as 



772(0;) = iv 



Py 



ai 



VpAx) 



^iS{xo,x)/h 



1 



a2 



\/Px{x)^jG{x) 



-iS(xQ ,x) / h 



To obtain ( 33 1 , we used the equalities 
v{x) +px{x) 



\Pv 



v{x) ~Px{x) _ ^ 

\Py\ " G' 



(33) 



(34) 



Inserting 771 (x) and 772 (x) into the current conservation law ( 26 1 , we find 

1 / 1 



\mixr-\mixW 



Px{x) \G{x) 
\Py\ 



G{x)\{\a,\'-\a2\'), 



\ai\ 



(35) 



where we assumed that the action 5(xo,x) is purely real. 

When we consider scattering from an electron region on the left to an electron 
region on the right, we can introduce the coefficients ai — 1 and a2 = r on 
the left, and ai = t and 02 = on the right. Since v — —\ on both sides, 
equation ([35|) tells us that 

(36) 



Therefore r and t can be treated as the reflection and transmission coefficient 
respectively. 

Now let us turn to scattering from an electron region on the left to a hole 
region on the right. From equation ([s]) one infers that holes with positive velocity 
Vx — OLq {px,x)/dpx have negative momentum, see also [H]. Therefore the 
coefficient of the right- moving hole state is 02, and we set a2 = t and ai = 0. 
Since v{x) > 0, we also have v = +1. Inserting this into the current conservation 
law (26), we again obtain equation (36 1. 

In a previous paper jl4j . the authors constructed the asymptotic scattering 
states with the help of a different method, and defined expansion coefficients 
that slightly differ from those given above. The relation between these two 
expansions is given in [Appendix C| 
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Figure 3: Stokes diagrams for the three different regimes outhned in section [T] a) Klein 
tunnehng, b) above-barrier scattering and c) conventional tunneling. Bold points show the 
turning points, the solid lines correspond to anti-Stokes lines and the wavy lines designate 
branch cuts of the function (z — zq)^!'^. In the figure u{z) = —z^. 



3. Stokes diagrams and the WKB approximation in the complex 
plane 



In the previous section we established that the functions -q^ (z), given by (21 ), 
solve the reduced equation ( 10 ) for 771 up to order h. Let us now choose the 
constant zq to be a turning point, i.e. a point where Pxi^z) vanishes; Pa;(^o) = 0- 
There are lines, called anti-Stokes lines [18, 19, 23J, emanating from the point zq 
along which the imaginary part of the function s(zo, z), given by equation (20), 
vanishes, that is Im[s(2;o, -2)] = Along an anti-Stokes line both asymptotic 
solutions are of order one with respect to the small parameter h. On the anti- 
Stokes line 7 the exact solution 'ip{z) can be represented as 



^P{z)^Clr,+ iz) + Clv^{z), 



(37) 



where the superscript '7' relates to the anti-Stokes line. Equality (37) only holds 
up to order h. 

Let us introduce a so-called Stokes diagram, in which the anti-Stokes lines 
are drawn in the complex plane, together with the choice of the branch cuts 
of the square root in the definition of Pxiz). In figure [sj we show the Stokes 
diagrams for the regimes outlined in section [T] The bold circles depict the 
turning points, the solid lines depict the anti-Stokes lines, and the wavy lines 
represent the branch cuts. Near a simple turning point, which means that Px{z) 
has a simple root z = zq, we can approximate Px{z) by a{z — zq)^^^, where a 
is some constant. It is then easy to show [TBI [TOl [53] that three anti-Stokes 
lines emanate from it. 

In figure [sja) , corresponding to the regime of Klein tunneling, one sees four 
real turning points. There are classically forbidden regions between the left two 



-'We remark that what we call an anti-Stokes line is sometimes denoted as Stokes line, see 



e.g. 
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and between the right two turning points, while the anti-Stokes hne connecting 
the middle two turning points represents the classically allowed hole region. 
The position of these turning points depends on the transversal momentum \py \ . 
When it goes to zero, the two turning points on the left (and on the right) come 
close together and eventually merge, that is, the classically forbidden region 
disappears, as discussed in [T3]. When \py\ becomes larger, the turning points 
on the left (and on the right) move further apart, and whenever uq ~ E = \py\, 
the middle two turning points merge. If \py\ increases further, they disappear 
off the real axis, and there are only two real turning points left. However, in 
this case Pxiz) acquires two complex roots, and hence we have two complex 
turning points. This is the situation in the conventional tunneling regime, the 
corresponding Stokes diagram is shown in figure [sjc) . We use this term because 
the situation is similar to that of a Schrodinger particle that tunnels through a 
potential hump. Obviously, the analogy can be used only if the complex turning 
points are sufficiently far from the real axis. 

In the regime of above-barrier scattering, figure [sjb) , all four turning points 
are complex. Since the potential u{x) is real on the real axis, the turning points 
come in complex conjugate pairs. Each of the two turning points closest to 
the real axis gives rise to one finite anti-Stokes line, and to two infinite anti- 
Stokes lines. If the potential u{x) vanishes along the real axis at |2;| — >■ oo, the 
infinite anti-Stokes lines approach horizontal asymptotes. If u{x) is unbounded 
at |x| — >■ oo, as in the figure |3]b), the infinite anti-Stokes lines approach the 
real axis. When \py \ becomes smaller, the upper (and lower) two turning points 
come close together and eventually merge whenever \py\ vanishes. When \py\ 
increases, the distance between points closest to the real axis becomes smaller, 
and when uq — E = \py\ they merge. When \py\ grows further, we once again 
end up in the conventional tunneling regime, figure |3]c). 

Now we can reformulate the scattering problem in terms of the Stokes dia- 
gram. Every scattering problem can be reduced to the determination of the 
coefficients C^^ , CZ^ along the anti-Stokes line 71 emanating from the turning 
point zi , under the assumption that the expansion coefficients , CZ along the 
anti-Stokes line 7 emanating from the turning point zq are known. Generally 
speaking zq and zi do not coincide. The problem of establishing the connection 
between the expansion coefRcients at different anti-Stokes lines is known as the 
connection problem. It was first found by Stokes [35], and was further elabor- 
ated by many others, see e.g. [HI There are various ways to solve it. The 
approach that was used to produce the first connection formulas [SHI 133 13H] : and 
is usually taken in textbooks on quantum mechanics, e.g. |39| . is to approximate 
the potential near the turning point, and to solve the resulting equation exactly. 
In its most rigorous form, this method is known as the method of comparison 



equations [281 ES [301 (El EB ISS] , and is summarized and applied in Appendix B 



In the remainder of this section we introduce a different approach, that was pi- 
oneered by Zwaan gS], and further developed in [2S [HI [^ [201 HH US] ■ In this 
method one passes from one anti-Stokes line to another along a suitable path 
in the complex plane avoiding the vicinities of turning points. In the rest of 
this section we give an introduction to this method, more details can be found 
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Figure 4: The Stokes diagram for a simple turning point zq. The wavy line depicts the cut. 
The blue arrows show the direction of the growth of the action s{zo,z). The letters 's' and 
'd' indicate the sectors where rj'^ is subdominant and dominant respectively 



in [Appendix A[ 

Following the idea first set forth by Furry [27], see also [l8 l 120 ] [T9 l [23] . we 
first consider the transition between two anti-Stokes lines emanating from the 
same turning point zq. We notice that though rj^iz) and % {z) are asymptotic 
solutions for (27), every single of them does not provide an asymptotic solution 
of any exact solution within an area around zq. This already becomes clear 
if we take into account that rjf have cuts emanating from zq, while the exact 
solution does not have such cuts. These cuts correspond to the definition of 
Px{z) in the complex plane. To resolve this apparent contradiction, let us first 
consider an exact solution r]i{z) which has an asymptotic expansion ri^{z) along 
an anti-Stokes line 7. We assume that the action s{zq,z) grows along 7 (as 
indicated by a blue arrow in figure [4| . For convenience let us choose the cut 
in the definition of Pxiz) along the anti-Stokes line next to 7 in the clockwise 
direction. If we now leave 7 in the counterclockwise direction, the action s(zo, z) 
acquires a positive imaginary part and ti^{z) gets exponentially small. Along 
the anti-Stokes line 71 we thus recover an asymptotic expansion in the form of a 
single "incoming" wave, meaning that s{zq, z) decays along this line. The 27r/3- 
sector between 7 and 71 is designated by 's' (subdominant) in figure |4j meaning 
that within it rf^{z) gets exponentially small values. In the next 27r/3-sector 
in the counterclockwise direction, designated by 'rf' (dominant), rj'^{z) becomes 
exponentially large. It still gives the asymptotic expansion of 771 (z) everywhere 
apart from some vicinity of the anti-Stokes line that coincides with the cut. In 
this vicinity the asymptotic representation fails since in the considered sector 
we can not avoid the appearance of the additional term given by r/^ (z) in the 
asymptotic of r?i(z) against the background of exponentially large values of 
r}i{z). Thus rj^{z) gives only one term of the asymptotic expansion of rii{z) 
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a) h) 



Figure 5: The Stokes diagram for two simple turning points zq a-nd ^i- The wavy Unes depict 
cuts. Blue arrows show the direction of the growth of the action s{zq, z). The letters in circles 
are given with respect to zq. Diagram a) corresponds to »7i(z) = rj'^ (z) along 7 and diagram 
b) to r?i(z) = T)~{z). 



along the right lip of the cut. Left and right lips of the cut are defined with 
respect to an observer standing on the cut with the turning point behind him. 
A similar consideration shows that the other term of the asymptotic expansion 
of rii{z) is given by 'q'l [z) on the left lip of the cut. 

For any point z on the cut we can define a nearby point on the right lip 
Zr and a nearby point on the left lip Zj. T hen one has Pxizi) = e~™ Px{zr)-, 
g^/'^{zi) — i'g~^/'^{zr), see also Appendix A Therefore we obtain the following 



asymptotic representation of the exact solution 771 (z^), i-e. on the right lip of 
the cut, 

m{zr) = 4^jMe-(-o,..)/'> + 5_^^(Me— (-o,..)A. (38) 

Px {Zr) Px {Zr) 

The latter can be written as 

m{Zr)^'ni{Zr)+W'ril{Zr)- (39) 

The fact that 771(2:) can not be approximated by r]'^{z) along every anti- 
Stokes line is associated with the Stokes phenomenon |35j . The constant i in 
front of 7yj~ {z) in p9| ) is called the Stokes constant. In contrast to the Schrodinger 
equation p]51ITOll20l[^ we found an extra factor v in front ri^{zr), which relates 



to the additional amplitude factor g{z), see Appendix A 

Now we turn to the case where the anti-Stokes lines emanate from differ- 
ent turning points. When we are dealing with a finite anti-Stokes line that 
connects the two turning points, the initial expansion is valid along the entire 
line. However, we do have to change the reference point of the action, that is, 
write s{zo,z) = s{zo, zi) + s{zi,z). This introduces the additional phase factor 
gis(zo,zi) -^yjjich the coefficients have to be multiplied. 

When we are dealing with two turning points zq and zi that are not con- 
nected by an anti-Stokes line, the situation is more complicated. In figure |5] 
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Figure 6: The potential u{x), togetlier with the effective potential Py — v^{x). There are two 
classically forbidden regions, separating the classically allowed electron and hole regions. 



one sees two simple turning points and corresponding anti-Stokes lines. Let us 
consider the same asymptotic solution 'r]^{z) as was considered above. In the 
region A4 bounded by four anti-Stokes lines, rif{z) gets dominant when viewed 
from the reference point zq. Against the backg roun d of exponentially large val- 
ues of ri^{z), an exponentially small term in (39), given by rj^iz), should be 
neglected within the accuracy of the WKB approximation. This means that 
the connection procedure prescribed by WKB-method is not bijective: func- 
tions differing by rj^ (z) on the clockwise lip of the cut will be mapped into the 
same expansions along the anti-Stokes lines emanating from zi. This leads to 
the so-called one-directional nature of the connection formulae, see [Appendix A| 

and [niiiniiMiin]- 

Let us summarize this section. We have formulated a set of rules on how to 
pass between different anti-Stokes lines. Using these rules one can solve certain 
kinds of scattering problems. However degenerate turning points cannot be 
treated within this approximation. In the next section we consider necessary 
generalizations of the WKB-method. 



4. Tunneling through a barrier supporting hole states 

In this section we solve the scattering problem for the first case of section [l] 
that is, tunneling through a barrier supporting hole states. The classically 
allowed region for this case has been extensively studied by the authors in [2] 
with the canonical operator method [M], and particular emphasis was placed 
there on its geometric interpretation. In this section we are mainly interested 
in the transition through the classically forbidden region. 

In figure [6] we show the potential u{x) of the original (Dirac) equation Q 
and the effective (classical) potential p'l — v'^{x) in equation (10). One sees 



that there is a classically allowed hole region, separated from the two classically 
allowed electron regions by two classically forbidden regions. Note that this 



13 



classically allowed hole region corresponds to the anti-Stokes line in the middle 
of figure |3^. We assume that this region is broad enough to use the semiclassical 
solutions ( |21[ ) within it. Therefore we can split the problem of transmission 
through a barrier supporting hole states into two simpler problems: 

1. Transmission from the electron region to the hole region; we will call this 
transmission through an n-p junction. 

2. Transmission from the hole region to the electron region; which will be 
denoted by transmission through a p-n junction. 

We start our treatment by introducing the transfer matrix, that connects 
the expansion coefficients 01^2 in the electron and hole regions, and relate its 
elements to the reflection and transmission coefficients for the n-p and p-n junc- 
tions. We proceed by obtaining the reflection and transmission coefficients from 
the complex WKB method. The formulas we find this way are not valid for 
near-normal incidence, that is, when the transversal momentum \py\ is small. 
Applying the comparison equation technique, [Appendix B[ we finally obtain 
expressions for the reflection and transmission coefficients that are uniformly 
valid in the entire range of incidence angles. 

4.I. Transfer matrix 

Let us start by considering an n-p junction, with an electron region on the 



left, and a hole region on the right. In these regions we have the expansions ( 30 ). 



We introduce coefflcients af, af and a'^, af corresponding to a right /left moving 
electron and right/left moving hole respectively. We then define the transfer 
matrix T„p as the matrix that connects the expansion coefficients (aj?,aj') and 



^-^^"pf i V ), (40) 



af j -""M a!' I' '"''"I T21 T22 

where the coefficients on the left are defined with respect to the turning point 
v{x-) = —\py\, and the coefficients on the right are defined with reference 
to the turning point a;+, where f = \py\. Now let us express the coefficients 
of this matrix in terms of the reflection and transmission coefficients. 

To determine T„p it is enough to know two linear independent solutions. As 
it was shown in section [2j one can take rji and 773. For rji we have 



r [ T21 T22 [ 



(41) 



Having already constructed r]2, equation (33), we obtain 

where the minus sign on the left comes from opposite values of v for the electron 
and hole regions. Since complex conjugation interchanges e*"^/'' and e^*"^/'', the 
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coefficients in (42) are swapped as compared to (41 1. Solving (41), (42), we find 
tliat 

1/t -r*/t* 

r/t -1/r 



T — 



(43) 



To find transfer matrix Tpn for a p-n junction, 



T 



T 



Til Ti2 
T21 T22 



(44) 



we exploit the fact that an n-p junction for rji is a p-n junction for r/2, a fact 



that follows from equation (10). Therefore two linear independent solutions in 



this case are 772 and r/l. Starting from the reflection and transmission problem 
for r]i , we find the following equation for transfer matrix Tp^ by considering 77* , 



1 



Til T12 \(t* 
T21 T22 M 



(45) 



In obtaining the above result, note that the right-moving hole wave is propor- 
tional to e"*'^/'', and that the electron and hole regions for rji are interchanged 
as compared to rji. On the other hand, from the solution given by 772 we obtain 



Til T12 
T21 T22 



From ( 45 ) , ( 46 1 we obtain the transfer matrix as 



T — 



1/t* -r/t 
r*/t* -l/t 



(46) 



(47) 



Formulas (43) and (47) show that once we have the reflection and trans- 



mission coefficients for an n-p junction, we immediately know the full transfer 
matrix for both an n-p and p-n junction. 

Now we come back to our initial problem, transmission through a barrier 
supporting hole states. The transfer matrix for this problem can be obtained 
from the two transfer matrices considered before, but one has to take into ac- 
count that the wave functions in the hole region are defined with respect to 
different reference points. Therefore we will need the additional matrix 



J. 7)' 











L ■ 



Pi 



(48) 



where 3^14. and X2- are the turning points on the left and on the right of the hole 
region, respectively. To distinguish the two different junctions we also denote 



T 

J 7 



pn 



i/i* -f/i 



Then we compute the total transfer matrix as 

T — T T T 

npn — -^np-^pp^pn- 



(49) 



(50) 
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As a; — >■ CXI, we have aj;' — t, af — and for x — > — oo we write = 1, 



af = r. This gives 



1 \ _ t 



I — ^npn 1 Q 



— Tnpn ■ (51) 



From (51) and the assumption that the transmission coefficient is real, cf. [22], 
we find the total transmission as 

1 _ tte-'-^/'' 
" 1 - ^*Pe-2'V'> ' ^^^^ 

which is the familiar Fabry-Perot formula [131 HI • 

^.2. Transmission coefficient from the complex WKB method 

Now let us analyze the problem of reflection and transmission using the 
theory outlined in section [3] We consider an n-p junction, and assume that 
the classically forbidden region is broad enough to allow for the use of the 
semiclassical wave functions between the two turning points, x^ < 

On the right-hand side, we start with the transmitted wave. 



7?i(x) ^ tJ '^^^\ exp ( - ^ / Px{x)dx\, (53) 



Px{x) \ h 
where Px{x) — ^v'^{x) ~ Py, and 

v{x) + ^/v^{x) — pi 

- — h — ^ 

We now choose the analytic continuation of the square root such that 

{v\x)-pI)'/^^^v^x)-pI x>x+, (55) 

so in this region G{x) = g{x), and rii{x) coincides with jyf (2), as defined in 
equation (23). Hence the analytic continuation of (53) reads 

Vii^)-t ^ /;^'(") expf-1 riv^z')-piy^'dz] (56) 



'(z;2(z)-p2)i/4 ^1 h 



Now we consider the behavior of rj^ along a certain path and compare it 
to that for the exact solution. The path is chosen by the requirement that the 
outgoing wave becomes subdominant when we move away from the real axis at 
positive infinity. Finally the branch cuts are chosen in such a way that the path 
does not cross them, see figure [Tj Following the green path, we first arrive on the 
anti-Stokes line where the exact solution is still accurately represented by 771(2). 
According to the arguments from section[3] this function then becomes dominant 
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a) b) c) 

Figure 7: Stokes diagrams for a) n-p junction, b) p-n junction and c) propagation through 
a single hump potential for a Schrodinger particle. Bold points, solid lines and wavy lines 
depict turning points, anti-Stokes lines and cuts respectively. The colored lines show the 
paths transferring the transmitted wave on the right of the classically forbidden region to the 
incoming (blue) and reflected (pink) waves on the left. The qualitative difference between the 
first two diagrams and the third one, i.e. between Dirac and Schrodinger particles, is that to 
transfer the outgoing wave to the incoming wave, one has to circumvent both turning points 
in the same direction in the Dirac case, whereas for the Schrodinger particle both turning 
points are circumvented in different directions. 

between the two turning points, and we can proceed with it. Considering the 
behavior of rjilz) on both Hps of the cut along the negative real axis, we can 
find the coefficients of the incoming and reflected waves. 

For an n-p junction, rii{z) becomes subdominant in the lower half-plane, so 
we choose the contour shown in figure [7]a). Then we have for x between x_ and 

{v^ - plf'^ = e""/yp2_„2(2.)^ x_ < x < x+. (57) 
Proceeding along the same contour, we find that on the lower lip of the cut 

{v^-plf'^^e-^-^v^{x)-pl, x<x^. (58) 

Therefore, the process of analytic continuation turns the outgoing hole wave 
into the incoming electron wave, as can also be seen from figure [7]a). One sees 
that both turning points are circumvented in the same direction. This is in stark 
contrast with the scattering problem for a conventional Schrodinger particle: in 
that case the same path transforms the outgoing wave into the reflected wave, 
and one has to circumvent both turning points in opposite directions to turn 
the outgoing wave into the incoming wave, see figure [7]c) . 

We just saw that upon analytic continuation along the green contour in 
figure [7]a, Pxix) turns into —px{x). However, we see from the definition of g{z) 
and equation ( |31[ ) that at the same time G{x) turns into —1/G{x), since v{x) 
changes sign. This means that the ratio G(x) /px{x) turns into 1/G{x)px{x), 
without an additional sign. To determine its phase on the lower lip of the cut, 
we can consider the limit py 0. Then in the lower complex half-plane we have 
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Figure 8: The Stokes diagram for an n-p junction. The contour C shows the path along which 
the integration should be performed in ||63|l. The dashed line shows the equivalent path. 



(v^iz) — Py)^^"^ — >■ v{z) and g{z)/px{z) '^/\Py\- Hence the process of analytic 
continuation gives the incoming wave on the left as 



'71(2^) 



t 



exp 



Px{x')dx' 



which can be rewritten as 



mix) 



texp{K/h) 



^Jpx{x)G{x] 



exp 




Px{x)dx 



where 



K = 



\Jpy — v^{x)dx. 



(59) 



(60) 



(61) 



Since the amplitude of the incoming wave should equal one, we find the trans- 
mission coefficient as 



t 



-K/h 



(62) 



Following the reasoning in 



] , we note that equality ( 62 1 holds for any 



distance between the two turning points. To show this, we note that instead of 
the contour shown in figure [7]a, one can take an arbitrarily large half-circle in 
the lower complex plane. As a result, the integration in equation (59) should be 



performed along the contour C shown in figure [Sj giving rise to the transmission 
coefficient 

t = exp (^^J^ J Px{z)dz^ , (63) 



Once we have obtained equation ( 63 ) , we can deform the contour in the complex 
plane to reduce this result to equations (62) and (61), see figure [s] 

In accordance with the theory from the previous section, we have to take 
the pink contour in figure [7]a) to reconstruct the reflected wave. Since we now 
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go around the left turning point in the counterclockwise direction, we find that 
on the upper lip of the cut 

iv'~piy^'^,Jv^x)-pl x<x^. (64) 

Studying the analytic continuation of G(x), we see that it turns into —G{x) on 
this lip of the cut. Therefore the reflected wave is given by 

771 (a;) = texp{K/h) (^-^J^^ cxp ^--^ J p^{x)dx^ . (65) 

Since we now have a negative number under the square root, we have to determ- 
ine its phase. To this end we consider how the incoming wave is transformed 
into the reflected wave along the circle circumventing the left turning point. 
Since the direction of this path is counterclockwise, Px{x) on the lower lip of 
the cut turns into e^'^p^ix) on the upper lip of the cut. Similarly, 1/G{x) is 
transformed into G{x) on the upper lip of the cut. Therefore the reflected wave 
becomes 



?7i(a;) = -iiexp(i4://i) W^^exp / p^{x)dx ] , (66) 



and with the help of equation ( 62 1 we conclude that the reflection coefficient 
equals 

r = -i. (67) 



The modulus of the reflection coefficient can be refined from ( 36 ) . We obtain 



|r| = - \t? = yr^^F^^, (68) 

which gives the uniform approximation for the modulus of the reflection coeffi- 
cient. Combining these two equations, r can be written as 

r = e-''^/2+«fl ^1 _ e-2K/h^ (69) 

where 9 is an additional phase factor. From our previous considerations, we 
conclude that 9 is small whenever both turning points are substantially separ- 
ated, i.e. when \py\ is not too small. The method we used until now fails to 
give a uniform approximation for the phase 9. This is related to the fact that 
any contour that transforms the outgoing wave into the reflected wave passes 
through the region between the turning points, if we do not want it to cross the 
cut. When both turning points become nearly degenerate, r]i(z) in this region 
ceases to be semiclassical. In the next subsection we will use the method of 
comparison equations to obtain a uniform approximation for the phase 9. 



With the results (62 1 and (69), we can write down the full transfer matrix 



for an n-p junction using equation (43 1. The determinant of the transfer matrix 
equals 

detT„p=^^. (70) 
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The approximation (67) gives detT„p — 0, so that T„p cannot be inverted. 



Therefore the transfer matrix obtained by using r — —i \s unidirectional. The 
unidirectionaUty closely relates to the fact that we neglected the subdominant 



solution between the two turning points. Using expression ( 69 ) for the reflection 
coefhcient, one obtains detT„p = — 1. Thus, the transfer matrix obtained with 
the use of equality ( 69 1 can be inverted. It is therefore bidirectional. 

Using equation (47), one immediately finds the transfer matrix for a p-n 



junction. Alternatively, one can derive the same result using the contour shown 
in figure [7]b). In Appendix A we show how to derive these matrices using a 



different formulation of the complex WKB method. 



4.3. An n-p-n junction and the method of comparison equations 

To obtain the transmission coefficient for a full n-p-n junction, we now 



use equation (52). Inserting the WKB result (67), we obtain a transmission 



coefficient that diverges at the transversal momenta satisfying the semiclassical 
quantization condition 



1 



^1+ 



dx' \J v'^{x') 



pI^tt \n + 



(71) 



where 3^14. and X2- are the left and right turning points at the border of the hole 
region, respectively. This divergence is due to the one-directional nature of the 
transfer matrix, as explained in the previous subsection. If we use ( 69 ) instead 



of (671 for the reflection coefficient and (62 1 for the transmission coefficient, we 



obtain 



Knp/h^ — Kpn/h^—iL/h 



-2K, 



ip/h -^y^ ^ — 2Kpn/h^ — 2iL / h-\-i7r—i0np — i6pri 



(72) 



where the quantity L is the classical action in the hole region, given by equa- 
and Kpn are the action integrals in the classically forbidden re- 



tion (48), Krip 



gions for the n-p and p-n junction respectively, and are given by equation (61 ). 
Finally, 9np and 9pn are the phases of the reflection coefficients (|69|) . For angular 



scattering, formula (72) gives a rather good result for the transmission coeffi- 
cient even if we put 6* = 0, as we show numerically in section [sj see also [H]. 
However for nearly normal incidence this result is no longer accurate. Thus the 
final step in the construction of the uniform approximation is to find phase 9 in 



(69) 



In order to obtain a uniform approximation for the reflection coefficient, we 
use the method of comparison equations. In Appendix B.3 we explain how to 
map the potential of a general n-p or p-n junction to a quadratic potential. The 
latter can be solved explicitly, and with the help of the mapping we can then 
construct an approximate solution of the original equation. We find that 9 is 
given by 



9 = Arg 



r 1 



K 
ttIi 



TT 

4 



K 

■nh 



K 
nh 



In 



K 
nh 



(73) 
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Using the asymptotic expansion of the F-function [JDJ SI] , one easily finds that 
9^0 when K/h is large, which agrees with the result of the previous section. 



Equation ( 73 ) was already anticipated in |14) , where it was obtained by explicitly 
solving the case of a linear potential, and then replacing the action between the 
two turning points by K. In|Appendix B.3]we now give a rigorous proof for this 



result. Further discussion of equation (72 1 is postponed until section [s] 

5. Semiclassical treatment of above-barrier transmission 

In this section we consider the second regime from section [T] namely above- 
barrier scattering. This regime can be split in two cases, i) scattering above a 
potential hump, and ii) scattering above a monotonous finite range potential. 

The first case describes for instance finite-range gating in graphene. The 
Stokes diagram for this potential is shown in figure [sjb). One sees that there 
are four turning points, two in the upper half-plane and two in the lower half- 
plane. In what follows we will assume that the potential tends to a constant 
at |x| — )■ oo. An example of such a potential is u{x) = exp(— x^). The Stokes 
diagram corresponding to this potential has infinitely many turning points. The 
approximation we made in section [3] is that only the four turning points closest 
to the real axis should be taken into account, while the others can be neglected. 

The second case is a monotonous finite range potential, that is, a finite 
increase of the potential. It corresponds to a single n-n junction, which can 
be used to model a transition between two macroscopically wide areas with 
different gates applied. The n-n junction can be simulated by the potential 
u{x) = tanh(a;). It is important to note that a finite increase cannot be captured 
by a finite polynomial, or a rational polynomial function. Therefore, the Stokes 
diagram corresponding to this potential differs substantially from that of case 
i), and should be considered separately. In section [7j we will construct the exact 
solution for this case, and also give a semiclassical treatment. 

Let us return to the Stokes diagram from figure [3]b) . In the first approxim- 
ation, we can neglect the two outermost turning points Z2±, that correspond to 
v{z±) = \py\, and consider only those in the middle, namely zi±, that corres- 
pond to v{zi±) = —\py\. The Stokes diagram we obtain in this way coincides 
with the one for a usual Schrodinger equation. The reflection coefficient for this 
case was first obtained in [421 143j , using a perturbation expansion. It was shown 
that the reflection coeffficient is exponentially small, in agreement with classical 
mechanics, where the above-barrier transmission is always equal to unity. We 
will use a different approach, that can be found in [HI |T51 |^ US]. In 
our derivation we will implicitly assume that the energy is comparable to the 
height of the potential hump and will not consider the transition to the Born 
approximation, i.e. to the case E ^ oo, which was studied in [44) . 

Let us first turn to the definition of the scattering states. When calculating 



the current (35), we assumed that the action s{zo,x) was purely real. However, 
in the case of above-barrier scattering the turning point zq is complex, which 
means that wc cannot take it as the lower limit of integration in the action, 
since it adds a complex part. Therefore, we take this lower limit to be a point 
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on the real axis. More specifically, we introduce so-called Stokes lines by the 
requirement that Re[s(zo,2;)] is zero, which implies that the action is purely 
imaginary. It turns out that the middle two turning points are connected by 
such a Stokes line, which therefore crosses the real axis, see figure [9] We call 
the point where this line crosses the real axis x^, and take it as the reference 
point for the action. 




C 



Figure 9: The upper half of the Stokes diagram for scattering above a short-range potential. 
The solid lines depict anti-Stokes lines, while the dashed line depicts a Stokes line, and the 
wavy line depicts a cut. The path in the complex plane is shown by the green and pink 
contour, and transforms the transmitted wave into the reflected wave. The integration contour 
C connects xq on the right lip of the cut with the same point on the loft lip of the cut and 
consists of the two black lines, the green line and the pink line. 

Let us now choose a convenient path in the complex plane. Starting at 
positive infinity, we first make a transition to the anti-Stokes line above or 
below the real axis. If the potential tends to a constant at plus (and minus) 
infinity, this transition does not change the expansion coefficients |45j , since the 
asymptotic solutions become exact and the distance between the real axis and 
the anti-Stokes line remains finite. 

To decide whether we should take a path in the upper or in the lower half 
plane, we start from the expected result. With respect to the reference point 
Xo, we have an incoming wave with unit amplitude, and we expect an exponen- 
tially small reflection. When changing the reference point to the upper turning 
point zi_|_, the coefficient in front of the incoming wave is multiplied by an ex- 
ponentially small factor, and the coefficient in front of the reflected wave is 
multiplied by an exponentially large factor, in such a way that they are of the 
same order of magnitude at the upper anti-Stokes line. This condition is neces- 
sary to pass between two anti-Stokes lines that emanate from the same turning 
point, see section [3j It turns out that this condition is not satisfied in the lower 
half plane. 

After these preliminaries, we can follow the procedure that was used for 
Schrodinger equations [SJ UHl \19\ [221 145] . We start with a transmitted wave at 
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positive infinity, 



where = ^Jv^{x) — and by Xq we denote the point XqzLe, where e oj^ 

We choose the analytic continuation of the square root as 

(V^X) - p2)l/2 ^ ^^2^^)_p2^ ^ > (75) 

{v^x)-piy/'^^^^v^x)-pl, x<xo. (76) 

which imphes that 771(2;) coincides with rj^{z) along the positive real axis. Let us 
consider ri^{z) along the path shown in figurejoj It becomes subdominant above 
the anti-Stokes line on the right. Therefore on the anti-Stokes line between 
zi+ and Z2+, the function rj^{z) correctly reproduces the behavior of 77(2), 
see section [3] When we continue along the contour, the transmitted wave is 
transformed into the reflected wave. Following the complex continuation, we 
see that Px{x) is transformed into e^'^px{x), and that 1/G{x) is transformed 
into G{x). Therefore, at negative infinity, we end up with 



X exp 



^v{x')^-pldx'\ , (77) 



Note that the integral in the second exponent is to be performed from Xq . The 
first integral is associated with the change of reference point; it goes along a 
contour C that connects x'^ with x^ . Since the square root has opposite signs 
on opposite lips of the cut, it can be rewritten as 

exp (^1 jy{z)-plY'^dz^ = exp (^^ £\v\z) - plf/^d^ . (78) 

Approximating the transmission coefficient by one, we find that the reflection 
coefficient is given by 

r = -^e^/^ K = 2i {v{zf ^ plf'^dz < Q. (79) 

The fact that if is a negative real number can be seen by performing the cal- 
culation for the prototype potential v"^ — p^ = z'^ + a^, see also 18, 19J. 



^Note that in the terminology of section [sj we would say that the point Xq is on the left 
lip of the cut, which is somewhat counterintuitive. 
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Figure 10: The upper half of the Stokes diagram for scattering above a short-range potential. 
The solid lines depict anti-Stokes lines, two of which are labelled by 71,2, while the dashed 
line depicts a Stokes line, and the wavy line depicts a cut. The green line shows the path 
taken in the complex plane. 



Approximation (79) does not hold when \py\ is close to E — uq, since in 



this case the middle two turning points are close together. We can get a more 
accurate prediction for the modulus of the reflection coefficient by considering 



the current conservation ( 36 ) 



In fact, equation (77 1 is equivalent to 



Combining this with |rp + |ip = 1, we find that 



f 



(80) 



(81) 



Note that we can now deform the contour into a large semi-circle, and as long as 
the upper-most turning point 2:2+ does not come into play, the above derivation 
still holds. Therefore we conclude that equation (81 1 holds regardless of the 



distance between the middle two turning points. The same result is obtained by 
applying the method of comparison equations, see [Appendix B.4 from which 
we also find the correct phases. 



J4' 



\/rTe2i?77^' 



Vl + e2^/^' 



where K was defined in equation ( 79 ) , and 



= Arg 





iK\ 


K 


K 






nh 


nh 



In 



nh 



(82) 



(83) 



All the above results coincide with the results for an ordinary Schrodinger equa- 
tion P] . 

Since the above results are identical to those for an ordinary Schrodinger 
equation, they fail to explain the total transmission at normal incidence. Indeed, 
when the transversal momentum \py\ becomes small, that is, for near-normal 
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incidence, the upper two turning points come close together (as do the lower 
two) and therefore they have to be treated as a cluster. As long as E — uq 
is sufficiently large, the two clusters can be treated separately. To derive the 
reflection coefficient, we now position the branch cut differently, and use the 
contour shown in figure [TOj 

To find the reflection coefficient in this approximation, we start with the 



outgoing wave (74) with t equal to one. Then we change the reference point, 



i.e. the lower limit of the integral, to zi^ and obtain 

= e^/2'>(z5)_^g..(.i+,.)/^ (84) 

Px (z) 

In |Appendix B.5| we use the method of comparison equations, to show that upon 
passing from the anti-Stokes line 71 to 72, the accurate representation of the 
exact solution becomes 

fj.iz) - e^/2'.(z5l_^^!Me-(-i+.-)/" _ ,„e/f/2/.(z|)^{!Me—(-i+. (85) 

Px (z) Px (z) 



where 



and 



r(-S'/7r/i) y -S/irh 



r^(S/7rh) \n{-S/Trh)~{~S/irh) 



(86) 



r-Z2 + 



S= / {v\z)-piy/^dz<0. (87) 
Changing the reference point back to xq, one sees that the first term in equa- 



tion (85) is the incoming wave normalized by one. The second term is the 



refiected wave, and the reflection coefficient equals 

r = -me^/''. (88) 

It is readily seen from the Stirling approximation for the F-function [IHl HI] 
that a approaches one when —S/nh becomes large. Upon normal incidence, the 
upper two turning points merge and S vanishes. One can check that a vanishes 



in this case, and therefore equation (88 1 correctly predicts total transmission 



at normal incidence. We see that similar to the case of conventional Klein 
tunneling, or scattering in the presence of hole states, total transmission comes 
from the merging of two turning points. However, this time they do not merge 
on the real axis, but in the complex plane. 



Equation ( 88 1 was derived under the assumption that the upper two turning 



points are close together, whereas equation (82 1 was derived under the assump- 
tion that the middle two turning points are close together. In the intermediate 
regime, we can combine the two expressions into 

\r\ = , , (89) 



25 



which shows the correct behavior for normal incidence. The transmission coef- 



ficient can be derived from this from the current conservation 



\t\' 



1. In 



section |8] we show that, surprisingly, equation (89) is in better agreement with 
our numerical result than equations ( 82 1 and ( 88 ) . 



6. Tunneling through a barrier without hole states 



Now let us consider the third regime from scction[TJ the conventional tunnel- 
ing regime. We consider a short-range potential, for which the Stokes diagram is 
shown in figure [2]c). Two of the four turning points are real, and the other two 
are imaginary. In the previous section, we saw that imaginary turning points 
give rise to exponentially small refiections, so we will start by neglecting their 
influence. 

To relate the transmission coefficient to the reflection coefficient, we use the 
contour shown in figure [7]c. Following the green and pink contours in a way 
similar to section [4j we find 



-K/h 



K 



dx 



v'^{x) > 0. 



(90) 



Note that in this case the analytical continuation of g{z), see equation (23) 



does not give rise to an additional phase factor. A second way to see this is 
by realizing that both turning points x_ and x+ correspond to v{x±) = — \Py\- 
Considering the additional amplitude factor in equation (21 ), and integrating it 
both turning points, we find 



exp 



v'{x)dx 
^pI - v{xY 



= 1. 



(91) 



Therefore, the result is essentially the same as that for an ordinary Schrodinger 
equation, and we call this regime the conventional tunneling regime. 



We note that equation ( 90 ) does not hold when the two turning points come 
close together. However, by similar reasoning as in section |4j see also [22] i we 
can conclude that the relation 



(92) 



does hold regardless of the distance between the two turning points, as long as 
there are no other turning points close to the contour in the complex plane. 



Combining equation (92 1 with the current conservation |r| + \t\ ~ 1, we can 



obtain the modulus of the transmission coefficient. In Appendix B.4 we outline 



how we can use the method of comparison equations to reconstruct the correct 
phases, with the result 



ci4> 



t = 



(93) 
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where (j) is defined by equation (83). Tliis result coincides with the one derived 
for above-barrier scattering, equation ( 82 1 , except for the fact that K is now pos- 
itive instead of negative. The same resuh is found for an ordinary Schrodinger 
equation [25] , 

In section [8] we compare our result to numerical calculations, and see that 
the agreement is reasonable. The discrepancy is due to the influence of the 
other two turning points, that were neglected in the above treatment. Finally, 
we note that equation ( 90 1 also holds for cases when there are more than two 
complex turning points, as long as these are not too close to the real axis, and 
the turning points on the real axis are not too close together. 



7. The exactly solvable model of the monotonous finite range poten- 
tial 

Up to now we have considered scattering of massless Dirac fermions by a 
potential hump. Implicitly, we assumed that we only have to keep four turning 
points in the Stokes diagram to find the major contribution to the scattering. 
This assumption naturally led to Stokes diagrams topologically equivalent to 
those for a parabolic potential, u(x) — —x^. In contrast to a potential hump, 
a monotonous finite range potential (finite increase) cannot be modeled by a 
polynomial function. This leads to a topologically different Stokes diagram. As 
an example let us consider an exactly solvable model for a finite range potential, 
provided by the function 

= Y(l + tanh(x)). (94) 

The exact solution for this potential was constructed in [33]. Here we present 
a slightly different approach, following the general method outlined in . 
Similar techniques was employed in [4 7) . where the eigenvalue problem for a 
potential well u[x) = — 1/ cosh(a;) was solved. 



Inserting the potential (94) into equation (10), we find the differential equa- 
tion ^ 

/i^^^ [g2 tanh^(a;) + qx tanh(x) + go] m = 0, (95) 



where 

92 = y (^y , qi^uoi^Y^^)' 9o = (^y - ^^j +«/jy • (96) 
To solve this equation, we first perform the substitution 

C = (1 -tanh(2;))/2, (97) 

leading to 

+ h'^ (1 - 20' + 91 (1 - 20 + go] m = 0. (98) 
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Substitution ( 97 1 maps the real axis to the interval < ^ < 1 in such a way the 
limit a; — >■ oo corresponds to ^ — >■ and x — oo corresponds to ^ — ^ 1. When 
^ — >■ the part in square brackets in equation (95) tends to pi, where 



Pi = ^ [uo - EY - p. 
and when ^ ^ 1, it becomes p|, where 



P2 = \jE'^ -pI- 
Therefore, let us make the substitution 



(99) 



(100) 



(101) 



After some algebraic calculations, one finds that w satisfies the hypergeometric 
differential equation. 



(1 - + (c - (« + + 1)0 ^-^bw = 0, 



with the parameters 

_ , »Pi m iuQ 



ipi ip2 lUo ipi 

b=— + — — — , c=l + — . 



2h 2h 



2h 



(102) 



(103) 



Two linearly independent solutions of (102) can be taken as formulas (15.5.3) 
and (15.5.4) in gD], 

wi=2Fi{a,b,c;0, = ^'"^(1 - ^)^"""''2Fi(l - a, 1 - 6, 2 - c; 0- (104) 
Hence rji can be written as 

r^i = ciri/2''(l - e)'^^/2^Fi(a, 6, c; 

+ C2r'^^^"'(l - e)"'P'/'Si^i(l - a, 1 - 5, 2 - c; 0- (105) 

When X ^ oo we can use the approximate relation ^ ~ e^^^. From the equality 
2Fi{a, 6, c; ^ = 0) = 1, we find the asymptotic behavior of rji 



rji — > cie ' + C2e ' . 
Therefore we conclude that the function 



(106) 



(107) 



where the superscript 't' stands for tunneling, gives the solution for the scat- 
tering problem in the regime of Klein tunneling {E + \py\ < uq). On the other 
hand, the function 



(108) 
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where the superscript 'a' stands for above, gives the solution for the scattering 
problem in the regime of above-barrier scattering {E — \py \ > uq). Now we use 
formula (15.3.6) from [40], 

F{a, b, c; = W^^^f^^^,2Fi (a, 6, a + 6 + 1 - c; 1 - ^ 
1 (c — a)i (c — o) 

r(c)r(a + 6-c)^^ 



mm 



(1 - <e)'""" 2i^i(c - a, c - 6, c + 1 - a - 6; 1 - 0, (109) 



that relates the values of the hypergeometric function at the singular points 
^ = and ^ = 1. Using that for x — )■ — oo, ^ ~ 1 — e^^, we find that 

(a) T{2- c)T{c~ a~b) .^^^/^^ V{2 - c)T{a + b - c) j^^^/h 
'1 r(l-a)r(l-fe) T{\+a-c)T{l + b-c) 

(t) r(c)r(e-Q-b) r(c)r(a + b-c) 

^ r(c-a)r(c-&)^ + r(a)r(6) ^ ' ^^^"^ 

To find the proper reflection and transmission coefficients, we remember that 
the scmiclassical scattering states are defined as (30). From these equations, 
we conclude that the ratio of the coefficients in front of Q-^^^x/h ^^^^ ^ik2x/h 
X — > — oo is equal to 



by I 
This gives 



(111) 



^(g)^ \Pv\ T{a + b~-c)T{l-a)T{l-b) 

E-p2T{l + a-c)T{l + b-c)T{c-a-by ^ ' 

\py\ T{a + b-c)T{c-a)T{c-b) 

E^P2 r{a)r{b)r{c-a-b) ■ ^ ' 

Analogous arguments give the transmission coefficients as 




e^Uq-p, r(i-a)r(i-6) 

E-P2 T(2- c)T{c - a-b)' 

\Pv\ r(c-a)r(c-6) 

P2 ^{Un-E+pi){E^p2) r(c)r(c - a - &) ■ 



(114) 
(115) 



Of course these transmission coefficients should equal one at normal incidence. 
To see this, one has to take the limit \py\ — 0, upon which pi — >■ \Uo — E\, 
P2 ^ E and E ~ p2 ~ Py/{2E). This means that for both the Klein tunneling 
regime and for above-barrier scattering, the factor in front of the F-functions 
becomes one upon normal incidence. Furthermore, we note that at py = 0, 
we have b = for the regime of Klein tunneling and c = 1 + b for above- 
barrier scattering. Therefore in both cases the quotient of F-functions equals 
one and the exact solution shows total transmission at normal incidence. By 
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A lm{z) 




Figure 11: Part of the Stokes diagram for scattering by a tangent hyperbolic potential. The 
cross at 2 = i7r/2 depicts the pole of the tangent, and 73 denotes the anti-Stokes line that 
starts at Z2+ and ends at the pole. Going along the contour C, and m aking use of the method 
of comparison equations, one finds the transmission coefficient Jllml. 



the same methods, one can construct the solution for the decreasing potential 
u{x) = uo[l — tanh(a;)]/2. 

To get more generic insight in the process of above-barrier scattering by 
a finite range monotonous potential, let us now treat scattering by a tangent 
hyperbolic potential semiclassically. In figure II one sees the corresponding 



Stokes diagram. The main peculiarity of this diagram is the existence of a pole 
at z = j7r/2, and the existence of a finite anti-Stokes line, that ends at the pole. 
To obtain an approximation for the transmission coefficient in the case when 
and Z2+ are far apart, we can apply the WKB approximation in the same way 
as in section [5j This takes only the right-most turning point into account, and 
we once again find the transmission coefficient (81). This answer however fails 
to explain total transmission at normal incidence. The reason for this is that 
at normal incidence: turning points zi+ and Z2+ merge, and the approximation 
considered above is no longer valid. 

To obtain a prediction that is also valid for near-normal incidence, one has to 
treat the two turning points in the upper half plane as a cluster. This problem 
can be solved by the method of comparison equations, see [Appendix B.6[ and 
gives the reflection coefficient as 



-ze 



where K is given by equation ( 79 1 , and 



'22 + 



(116) 



(117) 
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In turn, 9 is given by 
9 = Arg 



r 1 



iS 



S ^ ^ 
nh Tih \ Tih 



TT 

4' 



(118) 



Note that this result does give total transmission at normal incidence, since S 
tends to zero when \py\ tends to zero. 

The same answer can also be derived in a simpler way, by making effective use 
of our previous results for the Klein tunneling regime. Indeed, let us consider the 
half axis z — x + m x > 0. Along this axis we have tanh(a; + i7r/2) = coth(x). 
When X — ^ 0, the potential u{x) is proportional to 1 + coth(x) — > oo, which 
means that we are dealing with a p-n junction along this line! Along this line, 
the equation for rii{x + iir /2) reads 



pl+ihV'{x)\ri,=Q 



(119) 



where V{x) — C/o(l + coth(a;))/2 — E. Sufficiently far from the pole at a; = 0, 



we can use the transfer matrix for a p-n junction, equation (47), to establish 



the connection between the wave function on the anti-Stokes line 71 , on the the 
right of the classically forbidden region, and on anti-Stokes line 73, to the left 
of the classically forbidden region. This gives the relation 



pn , 1 1 



\/Px{x) 



-iS{x2^^x)/h 



+ Tp.n,2l 



^iS(x.2+,x)/h 



^iS(xi^,x)/h 



\/Vx{x)^jG{x) \/px{x)yjG{x) ' 



(120) 



where all quantities relate to equation (119), and xi^2+ = ^1,2+ ~ «7r/2. One 
can now establish the relation between G{x) in equation ( 120 1 and g{x + i7r/2), 
to find ?7i(2) in the complex plane. Since the term with coefficient Tp„ 21 turns 
out to be dominant between 73 and 72, it provides one of the two terms of 
the asymptotic expansion along 72. When continues to the real axis, this term 
becomes the reflected wave. Changing the reference point of the action from Z2+ 



to xo, one finally arrives at the previous result ( 116 ) for the refiection coefficient. 

Upon normal incidence, the turning points zi+ and Z2+ merge. Along the 
line X -\- m /2 we have conventional Klein tunneling. Therefore Tp„ 21 vanishes, 
i.e. there is no refiected wave. From the previous arguments, we then conclude 
that the reflection coefficient for above-barrier scattering vanishes as well. Thus 
total transmission for a particle that is normally incident on an monotonous 
increasing potential is related to conventional Klein tunneling in the complex 
plane. This effect can therefore be referred to as "virtual Klein tunneling" . 

Of course the same calculations can be done for a finite decrease of the 
potential. In this case one finds that the refiection coefficient is given by 



(121) 



where S* > is the action between the two complex turning points in the upper 



half of the complex plane, similar to equation ( 117 1. Given this definition of S 



the phase 9 is defined by equation (118) 
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Glueing the increasing and decreasing potentials together, we obtain a finite 
step potential of an arbitrary width. In contrast to the potential hump, con- 
sidered in section[5j the transmission through this structure reveals Fabry-Perot 
oscillations, 

\t I - l^ill^^l (122) 



where ri and r2 are given by equation (121 1 with parameters corresponding to 



the left and right junction respectively. The transmission coefficients are given 
by 1^1,2! = \/l — |?'i,2p- Finally, L is the action between the two reference points 
for the separate junctions, xo_ and Xo+, 



L = 



^v^{x)-pldx. (123) 



8. Comparison with numerical results 

In this section we compare our semiclassical predictions from the previous 
sections to numerical results. These are obtained by approximating the potential 
by a series of small steps. Since the potential is constant between each of 
them, one can use the exact solution for a constant potential [TU]. Matching 
the coeflficients at each interface with the help of a computer, we obtain the 
reflection and transmission coefficients. 

Let us start by considering a finite increase of the potential, which corres- 
ponds to an n-p junction for E < uq. We model it by 

u{x/li) = 0.5 Wo [1 + tanh(10x//i - 5)] . (124) 

When X changes from to /i, the potential saturates with an accuracy of 0.01%. 
Therefore, the junction can be cut at these points, without any substantial 
numerical error. An n-p-n junction is modeled as an n-p junction with length Zi, 
a p-n junction with length ^3 and a constant potential of length I2 in between. 



Figure 12 shows the transmission \t\'^ depending on the angle of incidence 0, 
which is related to the transversal momentum py by Py = Esin(f). We compare 
the numerical transmission for an asymmetric n-p-n junction with the semiclas- 



sical result, equation (72) with = 0, and the uniform result, where 9 is given 
by equation ( 73 ) . In the computation of the dimensionless parameters we take 
the maximum of the potential Uq and the length li/2 as the typical values vpo 
and I introduced in section [T] Then figure [12] corresponds to h = 0.08 and 
E = 0.4. In the case of graphene, this would correspond to an electron energy 
of 100 meV, a barrier height of 250 meV, and length scales = 70 nm, I2 = 150 
nm and ^3 = 90 nm. 



We see that the agreement between the numerical result and equation ( 72 1 



with 9 — becomes better as the angle of incidence increases, that is, deep in 



the semiclassical regime. Indeed, when we use equation (73) for 9, we uniformly 
approximate the numerical data over the entire range of incidence angles. Con- 
cerning the validity of the semiclassical approximation, we note that the agree- 
ment improves when the potential is smoother, i.e. when li and ^3 are large. 



32 



90 










<t>y 







15° 





Figure 12: Left panel: The transmission for an electron incident on an n-p-n junction as a 
function of the angle of incidence <p. The dimensionless parameters are h = 0.08 and E = 0.4. 
The barrier width I2/I = 4.3 and the n-p and p-n regions have lengths li/l = 2 and I3/I = 2.6, 
respectively. The blue line (solid) shows the numerical results for 99 steps, the green line 
(large dashes) shows the WKB approximation, equation \72\ where 9 = 0, and the red line 
(small dashes) shows the uniform approximation, where 9 is given by equation ( |73[ l. Right 
panel: Top view of the potential, showing the angle of incidence 0, and side view, showing 
the length scales ^1-/3. 



The exact solution obtained in section [7] perfectly coincides with the numerical 
results. 

One sees that apart from total transmission at normal incidence, there are 
also additional side resonances. This is a well-known phenomenon for transmis- 
sion through a metastable hole state. For a more detailed consideration, it is 



convenient to rewrite equation ( 72 1 in the form 



1 _ ^(1 _ e-2if.p/?»)(i _ e-2^fp^/^) + 4^(1 - e-2^"p/'')(l - e-^^p^/f") sin^ i) 

(125) 

where 'd = L/h+9np/'2+0pn/2 — TT/2. From this expression one immediately sees 
that the transmission coefficient is maximal at i? = nir, and that its modulus 
equals 

g-(A:„p+Xp„)//i 

'^"^"''■^^ = l_y(l_e-2^"./'')(l-e-2^fp./^ ^^^^^ 

For a perfectly symmetric junction, Knp — Kp^, the amplitude is unity, and 
the corresponding angles were called 'magic' (TUIIIII. For a generic asymmet- 
ric junction the height of the resonances decays. When we consider the truly 
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Figure 13: The angular depe ndence of the transmission for an electron incident on an short- 
range potential profile (128 i. The dimensionless parameters are h = 0.2 and E = 2. The 
length scale in the potential is given by li/l = 2. In both panels the blue line (solid) shows 
the numerical result with 49 steps, and the red line (small dashes) shows the semiclassical 
prediction for the transmission in the regime of above-barrier scattering derived from equa- 
tion l |88| l. In the left panel, the green li ne ( large dashes) shows the semiclassical result l |82[ l, 
that smoothly goes over in the result ( |93[ l for the conventional tunneling regime. In the 
right panel, the black line (large dashes) shows the transmission derived from the reflection 
coefficient (|89| for the regime of above-barrier scattering. 



To test our semiclassical results for scattering above a short-range potential, 
we use the model potential 



u(x/li) 



Uq 



cosh(10x/Zi - 10)' 



(128) 



Without substantial numerical error, we can cut this junction at a; = and 
X = 2li. In figure 13 we show the transmission for an electron incident on such 



a short-range potential. The dimensionless parameters are given hy h — 0.2 and 
E — 2. For the case of graphene, this corresponds to a particle of energy 200 
meV, a potential height of 100 meV and a length li = 70 nm. 

In the left panel of figure |13[ we compare the numerical transmission with 
the semiclassical result derived from (88) and the equality |ip = 1 — |rp, and 



the semiclassical result (82). Note that we use cartesian plots instead of an- 



gular plots from now on, in order to make the difference between the different 
approximations more pronounced. As anticipated in section [5j the semiclassical 
prediction (88), that takes the upper two turning points into account, works well 



at normal incidence. However, once we get closer to the point where the above- 
barrier scattering regime turns into the conventional tunneling regime (around 
30 degrees in the figure), the discrepancy becomes larger. Equation (82), which 
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was derived by considering the middle two turning points as a cluster, gives a 
slightly better result at this point. Note that this result for the transmission 



smoothly goes over in (93), which is seen to give a reasonable prediction for the 



conventional tunneling regime. The discrepancy is due to the influence of the 
two complex turning points, which were not taken into account in the derivation. 

In the right panel of figure |13[ we compare the numerical result with the 
results obtained from the semiclassical prediction (88) and the equality |ip = 
1— |rp, and our heuristic formula (89). We see that for above-barrier scattering. 



the heuristic expression ( 89 ) is more accurate than equation ( 88 ) . Overall, the 



agreement improves when the potential gets smoother, that is, when li increases. 
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Figure 14: The angular dependence of the transmission for an electron. In the left panel we 
consider a finite range increase of the potential, with h = 0.3, E = 2 and a length li/l = 2. 
The blue (solid) line shows the numerical result with 49 steps, and the red (dashed) line shows 
the semiclasical result ( |116| . In the right panel, we consider a broad potential barrier with 
h = 0.2 and E = 2. The width of the constant part I2/I = 2.9 and the potential increase 
and decrease have lengths li/l = 2 and (3 = 2.6, respectively. The blue (solid) line shows 
the numer ical r esult with 99 steps and the red (dashed) line shows the semiclassical result, 
equations (|122| and (|121|. 



To test our results for scattering above a finite increase of the potential, 
we go back to the potential (124). In the left panel of figure 14 we compare 
our semiclassic al res ult (116) for this case, with the numerical result and the 
exact solution (114). The dimensionless parameters are h = 0.3, E = 2 and 
li/l = 2. In the case of graphene, this corresponds to an energy of 200 meV, a 
potential of height 100 meV and a length li = 50 nm. The exact solution and 
the numerical transmission coincide exactly, and therefore only the numerical 
result is shown. We see that there is good agreement between the numerical 
result and the semiclassical prediction. Once again, this agreement improves as 
the potential gets smoother. 

In the right panel, we consider above-barrier scattering for a broad poten- 
tial hump that is modeled as a potential increase of length li, a constant part 
of length I2, and a decrease of length ^3. The dimensionless parameters are 
h = 0.2, E — 2 and li/l — 2. In the case of graphene, this corresponds to an 
energy of 200 meV and a potential of height 100 meV. The lengths are li = 70 
nm, I2 = 100 nm and Z3 = 90 nm. The semiclassical result is given by equa- 
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tion ( 122 1, with reflection coefiicieiit ( 121 1, and the modulus of the transmission 
coefficient is constructed from |rp — 1— Since the exact solution once again 
coincides with the numerical result, it is not shown. The agreement between the 
semiclassical prediction and the numerical result is quite good, and improves as 
the potential gets smoother. However, the positions of the maxima in the os- 
cillations are not perfectly reproduced. This is due to the fact that the turning 
points above and below the real line are quite close in this case, that is, Kjh is 
of order one, and our approximation is not strictly valid. 



9. Conclusion 

In this paper we have studied potential scattering of massless Dirac fermions 
in semiclassical approximation. We have shown that, depending on the energy 
of the incoming particle and its angle of incidence, there are three different re- 
gimes. These are i) the regime of the Klein tunneling, i.e. the regime when 
the scattering of electrons is mediated by hole states supported by the barrier, 
ii) the above-barrier scattering regime, and in) the conventional tunneling re- 
gime. For each of these regimes we found easy-to-use analytic expressions for 
the transmission and reflection coefficients. We have shown that the conven- 
tional WKB method does not allow to study near-normal incidence, due to the 
degeneracy of turning points at normal incidence. Therefore we cannot ob- 
tain expressions for reflection and transmission coefficients uniformly valid for 
any incidence angle. For near-degenerate turning points, the initial problem 
has to be reduced to a certain comparison equation with a well-known analytic 
solution. Using the solution of this comparison equation, we were able to ob- 
tain reflection and transmission coefficients for near-normal incidence, which is 
clearly crucial for physical applications. We completed the analytical part of 
the paper with the consideration of an exactly solvable model that simulates 
a monotonous n-n junction. This case is somewhat peculiar due to the pole 
of the potential in the complex plane. The unconventional Stokes diagram for 
this case is closely related to the one for the Klein tunneling regime. Therefore 
above-barrier scattering for such a junction can be treated as "virtual Klein 
tunneling" in the complex plane. The predictions provided by our analytic 
expressions show good agreement with numerical calculations. 
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Appendix A. The complex WKB method 



In this appendix we summarize the WKB approximation in the complex 
plane, as it was developed in[5Sl[2ni|IHl[2nilIlll3[lI]- We start with the 
explanation of the general method, and then discuss its application to potential 
scattering for massless Dirac fermions. 

Appendix A. 1. General formulation 
We start with the equation 

h^^+q{zMz)=0, (A.l) 

where /i ^ 1 is a small parameter and q{z) is an analytic function of the complex 
variable z, that may also depend on h. It has two approximate solutions, 

/i(zo, z) - exp (^1 £ dz' , 

Mzo, z) = exp (^-1 £ dz' q^l\z') 



(A.2) 



that will be referred to as basis functions from now on. Similar to the main 
text, we introduce anti-Stokes lines by the condition that 

s(zo,z) = r q^^^{z')dz' (A.3) 



is a real function. Here zq is a turning point, defined by the requirement that 
g(zo) = 0. On each anti-Stokes line 7 the exact solution can then be represented 
as 

4'iz)^CUi{zo,z) + C]f2izo,z). (A.4) 

The main problem when approximating the exact solution in this way is given by 
the Stokes phenomenon [35, 18l[Tnilini[13]i that we briefly touched in section [3} 



the exact solution V'C-z) has different representations (A.4) in different sectors of 
the complex plane. This naturally leads to the connection problem [18l [191 EHl 
[23] : given certain constants Cj, on the anti-Stokes line 7, which constants 
C7\ C]^ are needed to represent the exact solution on the anti-Stokes line 71? 
To connect the coefficients along the anti-Stokes lines 7 and 71, we introduce 
the matrix M, 

In the rest of this appendix we will determine the matrix M for the various 
transitions described in section [3] We will not give a precise estimate of the 
errors that are involved. Instead, we mention that precise estimates for the 
error terms are derived in [121 [13] . 

We start by considering the transition between two anti-Stokes lines that em- 
anate from the same turning point. Let us consider the situation of figure [A. 15[ 



, = M ^ ) . (A.5) 
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Figure A. 15: From a simple turning point (large circle) three anti-Stokes lines (solid) emanate. 
The wavy line indicates the branch cut, and the blue arrows indicate the direction in which 
the quantity s{zo,z) increases. The points zi and Z4 lie on the left and the right lip of the 
cut, respectively. Both 71 and 74 indicate the same anti-Stokes line, and are only used to 
distinguish different lips of the cut. The points 22 and 23 lie on the other two anti-Stokes lines 
72 and 73 respectively. 



with a simple turning point and a branch cut along the positive axis. Remem- 
ber that in section [3] we defined the left and right lip of the cut with respect to 
an observer standing on the cut with the turning point behind him. When we 
are on the right lip of the branch cut, we denote the anti-Stokes linealong the 
positive X-axis by 71. On the left lip of the branch cut, we use the notation 74. 
When proceeding from 71 in the clockwise direction, we first arrive at 72 and 
then at 73. Now assume that the expansion coefficients on 71 are given, and 
that s{z(j, z) increases along the anti-Stokes line, as indicated by the blue arrow 
in figure |A.15| This implies that when we move away from 71 in the clockwise 
direction, the action s{zq, z) obtains a negative complex part, as can be derived 
from the Cauchy-Riemann relations [18]. Therefore the basis function fi{zo,z) 
attains exponentially large values {dominant term), whereas f2{zo,z) becomes 
exponentially small {subdominant term). At a certain distance from 71, the 
subdominant term will be much smaller than the error in the dominant term, 
and we cannot keep the subdominant term within the accuracy of the method. 
Closer to the anti-Stokes line 72, both terms become comparable again. How- 
ever, the information about the coefficient in front of the subdominant term has 
been lost. The coefficient in front of the dominant term does not change, so we 
have = C^^ . In the most general form, the relation between the constants 
Cj^2 ^^'^ reads 

Let us consider two linearly independent solutions V'C-^) and 'ip{z), with coeffi- 
cients C7 2 and C7 2 on the anti-Stokes line 71 , respectively. This can be written 
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as 

(A.7) 

Taking the determinant on both sides, we see that the left-hand side is just the 
Wronskian, which is constant due to the current conservation for second order 
ordinary differential equations. On the right-hand side, the first determinant is 



also constant, which can be verified by inserting the definitions (A.2|. Hence 
the determinant of the second matrix is constant, and this constant does not 
depend on the anti-Stokes line 71 . We can then consider the transition from the 
anti-Stokes line 71 to 72, and write 



Cf Cf \ _ ( I 



^2 ^2 




{A.. 



Taking the determinant on both sides and using the fact that the determinants 
of the matrices with coefficients are equal, we find that the determinant of the 
first matrix on the left should equal one, and therefore /3 — 1. 

The change in the subdominant coefficient is therefore given by the so-called 
Stokes constant times the dominant coefficient, see also [151 HH 1201 [13j : that is, 



CJ' \ _ f I \ f 



-171 



(A.9) 



where a is the Stokes constant. The fact that only the coefficient in front of the 
subdominant term changes was called the 'principle of exponential dominance' 
in |21) . Note that when we start with a subdominant term only, its coefficient 



is unchanged, as can also be seen from equation (A.9). 

Now let us compute the actual value of the Stokes constant. Following [57], 
we start by noting that the exact solution should be single-valued when one 



makes a full turn around the turning point. However, the basis functions (A. 2 1 
contain the square root of z, which has a branch cut in the complex plane. Let 
Z4 be a point on 74, on the left lip of the cut, and let zi be the same point, but 



this time on the right lip of the cut, on 71, see figure A. 15 When we assume 
that we are dealing with a simple turning point, so that in the vicinity of the 
turning point, q{z) = a{z — zq) = re*"^, we can write 

q{z,)^re'', g(zi) = re^*+2", q^/'{zi) ^ e'^q'/^z^), (A.IO) 

where 6 is the angle at which the branch cut emanates from the turning point 
and equals zero in figure |^T5j We find that 

fi{zo,zi) = -if2{zo,Z4), f2{zo,zi) ^ -iJi{zQ,Zi). (A.fl) 

Since the exact solution is single- valued, 

V'(Z) =C7VlU,^l) + C2^V2(^0,2l) =C7VlU,Z4)+C2^V2(^0,^4). (A.12) 
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From equation ( A.ll I, wc then find that 



C7* = -iC]\ = -iCj\ (A.13) 

Hence, the Stokes constants have to be chosen in such away that when we go 
from zi to Z4 in the clockwise direction, see figure |A.15[ the matrix M from 



equation (A. 5 1 reads 



M^(_0^ -;). (A.14) 

From the above logic we know that fi{zo, z) is dominant in the region between 
71 and 72, and therefore the matrix A that connects the coefficients can be 



found from equation (A.9) 



Since /i(zo,^) is subdominant between 72 and 73, the matrix B that connects 
the coefficients on these anti-Stokes lines, is the transpose of the matrix in 



equation (A.9). Between 73 and 74, fi{zQ,z) is once again dominant, and we 



find the matrix C, 

^ 1 y ' V 1 

From the identity M = ABC, one then finds that 



B=(l ?V C^( I » ). (A.16) 



1 + &C b \ f -i 

a{l + bc) + c 1 + ah ) ^ \ -i 



(A.17) 



One of these equations turns out to be redundant, and solving the remaining 
three we find that all three Stokes constants are equal; 

a^b = c = -i. (A.18) 

When one goes in the counterclockwise direction, it turns out that all Stokes 



constants equal i. Since our basis functions (A. 2 1 are only accurate up to order 
h, we emphasize that the Stokes constant we derived here has the same accuracy. 
For a very precise error estimate, see [121 [23] ■ Finally, we mention that for a 
cluster of two turning points, from which four anti-Stokes lines emerge, a similar 
procedure one only gives a relation between the different Stokes constants, but 
not their actual value [25] . 

The second case, when two turning points are connected by a finite anti- 
Stokes line, was already completely discussed in the main text. We therefore 
turn to the third case, when we have a transition between two anti-Stokes lines 
7i and 72, that emanate from different turning points zq and zi. As discussed 
in the main text, one of the basis functions, say fi{zQ,z), for definiteness, will 
be dominant, and the other one, say f2{zo,z) will be subdominant. Under 



the assumption that and are comparable, it is shown in |23j that 
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sufficiently far from the turning point zq, the solution is accurately represented 
by the dominant term only, and the subdominant term should be neglected. 
Furthermore, an estimate of the error is derived. In terms of the matrices, this 
statement can be cast in the form 

/ is{za,zi)/h n \ f^l 

M=i^^ g JJj, s{zo,z,)^J dz'q'/^z'). (A.19) 

The above reasoning leads to the simple statement that [23] "one cannot proceed 
with an approximate solution, or an exact solution with approximately known 
initial conditions, in a classically forbidden region from the initial point in the 
direction in which the wave function decreases." This leads to the so-called 
'one-directional nature of the connection formulae' [TSJ [121 HSJ [H] ; a connection 
formula between a classically forbidden region and a classically allowed region 



can only be used in one direction. In equation (A.19) the one-directionality is 
manifest, since the matrix M has zero determinant. 

One may ask what happens when one does keep the subdominant solution 
between 71 and 72, assuming that the coefficient in front of it does not change. 
In [ini [13] ) the result of such a procedure was compared with the exact solu- 
tion for a parabolic potential, and it was shown that such a naive procedure 
gives wrong results for the exponentially small corrections. We therefore stress 
that the WKB-method applied to simple turning points can only give results 
in the leading-order approximation. When exponentially small corrections are 
required, one needs to resort to either unitarity arguments or one needs to make 
use of an exact solution, in the way that is explained in the main text and 
in Appendix B 



Appendix A. 2. Application to potential scattering for massless Dirac fermions 
We now want to make the connection between the abtract theory from the 
previous subsection, and the particular case of potential scattering for massless 



Dirac fermions which is considered in the main text. From equation (10 1, we 
see that for the case of graphene 

q{z)^v'^{z)-pl + thv'{z). (A.20) 



This makes the basis functions ( A.2[ ) slightly different from the asymptotic solu- 



tions (21) introduced in the main text. However, we can recover the latter by 



expanding the square root of q{z), 



{v^{z)-pl+^hv\z)f'' = {Az)-PIY" + ,(^2(y;l)V2 - (A.21) 



Since we only considered the leading order when constructing the functions ( A.2 ), 
corrections of order do not play a role. Therefore we can also neglect the 
correction in the amplitude factor and we recover the approximate solu- 



tions (21). Also note that the 'quantum' part in q{z) also slightly shifts the 



turning points, as compared to the way they were introduced in the main text. 
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However, one can show that this shift is of order h and that it also does not 
change the results to leading order. 

At this point we have to choose how we look at the term 




in the solutions (21). The first option is to look at this term as being part of 
the action. In that case, we have to change the lower limit in the integral when 
making a transition from one turning point to another. We see from the previous 
subsection that the Stokes constant for a clockwise rotation equals —i in this 
case. Note that to leading order this result is not changed by the correction 



in equation (A.20). Using the method of comparison equations, we give an 
independent proof of this fact in [Appendix B.2[ This first option naturally 
arises for the method of comparison equations, see [Appendix B[ 

The second option is to calculate the integral once, and to regard it as an 



amplitude factor. This was done in the main text, and leads to the solutions ( 23 1 



that include the additional amplitude factor g{z). However, in this case equa- 



tion (A. 11 1 is no longer valid. When v{zo) is positive, i.e. we are dealing 
with a turning point that separates a hole region from a classically forbidden 
region, one still has g^^'^{zi) = g^^^^{zi). However, when v{zq) is negative, 
i.e. we are dealing with a turning point that limits an electron region, one has 
g^/'^[zi) = — g~^/^(z4). These two statements can be combined as 

fi{zo,zi) ^ ~wf2{za,Zi), f2{zo,zi) ^ ~ii^fi{zo,Z4), (A. 23) 

where v — sgnw(a;o). Repeating the derivation presented in the previous subsec- 
tion, one finds that the Stokes constant for a clockwise rotation equals —iv, cf. 



equation (39). To comply with the main text, we will choose this second option 
in this appendix. 

We can now rederive the results from section |4.2| using the matrix approach 
that we just explained. The Stokes diagram for an n-p junction is shown once 
more in figure |A.16[ In section |4.2[ we had to carefully select the path that 
we take in the complex plane, but now we no longer have to be so careful. We 
therefore start on the anti-Stokes line 71 with two coefficients, Cj^ and C2^. 
Taking the path indicated in the figure, the analytic continuation of the square 
root is defined as 



iv'ix)- pt^y^' =^Jv^ix)-pl, x>x+ 

(w2(:E)-p2)i/2 =e-''^/^^pl-v^{x), x^<x<x+, (A.24) 



X < x_. 



With this definition /i(a;+, z) is dominant between 71 and 72, and since v{x-\-) > 0, 
we use the matrix ( A.9[ ) with a = — i, that is. 



-171 



(A.25) 
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Figure A. 16: The Stokes diagram for an n-p junction, together with the path we take in the 
complex plane. The solid circles represent the turning points x±, with v(x±) = it|py|, and the 
wavy lines represent the cuts. The relevant anti-Stokcs lines (solid lines) are labelled 71-4. 



Upon going from the anti-Stokes line 72 to 73, we see that i]^ is subdominant, 
whereas rj^ is dominant. Therefore we use the matrix 



73 







(A.26) 



In the region between 73 and 74 the function r\-^ {z) is dominant and since 
v(a;_) < 0, our matrix becomes 



1 



^73 
^73 



(A.27) 



The final step is to match the coefficients of the semiclassical solutions ( 23 ) to 
the coefficients afi and ^ of the right- and left-moving electron and hole waves, 
see section [i] Using that'c7^ = af,C$'^ = C7' = -af and Cj" = -af and 
multiplying the three matrices above, we find that for an n-p junction 



(A.28) 



where K is given by equation (61 1 . The result for a p-n junction can be derived 



from this by using the relation between the transfer matrices (43) and (47) 



The refiection and transmission coefficients that can be obtained from equa- 



tion (A.28 1, coincide with equations (62 1 and (67). We once again note that 



this matrix has zero determinant. This is a direct consequence of the fact that 
we neglected the exponentially small solution within the barrier, and implies 
that we can only use these matrices in one direction, "from the right to the 
left." 
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Appendix B. The method of comparison equations 



In this appendix we consider the second way to solve the connection problem, 
the method of comparison equations [2H1 [Ml [SHI El 132] ■ The basic idea of this 
method is to express the solutions of the original differential equation in terms of 
the solutions of a 'related' equation that is exactly solvable. This is the rigorous 
formulation of the approach usually taken in text books on quantum mechanics, 
see e.g. 1391 . where one approximates the potential by a linear function, and 
writes down its solution in terms of Airy functions. Generally speaking, the Airy 
equation is the prototype equation if a simple turning point is taken into account. 
If more turning points are considered, one can reduce the initial problem to more 
complicated prototype equations with the same number of turning points. In 
this appendix we use the formulation given in |32j . 

Appendix B.l. Explanation of the method 

We consider the second order differential equation 



,d^ 

dz 



2 + R{z,h)iP{z) =0, (B.l) 



where ft. ^ 1 is a small parameter, and z lies in a (possibly complex) domain D. 
The function R{z, ft) is supposed to be analytic, with asymptotic expansion 

oo 

i?(z,ft) = ^i?„(0)ft" (B.2) 

with respect to ft. A point zj at which Ro(z) has a root of order mj is called a 
turning point of order "^jj^ The total number of turning points is denoted by 
N + 1 and we set /i = J2j=o ™i- 



We will reduce (B.l ) to the related equation 



h^^ + Qicj,,h)V{cj,)^0. (B.3) 



The exact choice of Q{(l>, ft) will be specified below. Following jS^, we write 

V'(z,ft) = (0'(z))-i/V(0(z)). (B.4) 



Substituting this into equation (B.l ), we find that it is satisfied if 



•^Note that this definition is different from the one employed in [Appendix A[ where we 
called a zero of q(z) = R{z, h) a turning point, but that it coincides with the definition in the 
main text. 
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We will solve this equation order by order in h, assuming that 



oo 



n=0 

Collecting all terms of order we find that 
Gathering terms of order h^, we obtain 



Differentiating (B.8I and substituting the result into (B.9|, we obtain 



dz' R^^^^ 



(B.6) 
(B.7) 

(B.8) 
(B.9) 

(B.fO) 



We assume that the mapping (t){z) is non-singular, i.e. does not vanish 
within D. According to equation ( |B.8 1, this means that Qo{<f>Q) should van- 
ish whenever Ro{z) vanishes. By differentiating equation (B.8), one sees that 
Qo{<po) should have a root (poi^Zj) of order ruj at every turning point Zj. Thus 
we conclude that Qo{4>o) and Ro{z) have the same number of turning points 
within D, and that the order of their degeneracy coincides. 

In the simplest cases, we can choose Qaif/)) to be a polynomial function. 
Since we just concluded that Qo('/'o) should have the same number of turning 
points within Z?, we can write this polynomial as 



Qo(0) 



N 

J=0 



(B.ll) 



Taking the square root of (B.8 1, substituting the above expression and integrat- 



ing from a turning point to an arbitrary point z, we find 

dz'l^jRoiz')]'/'. 



dsl[\ 

00 (ZO) j^Q 



(B.12) 



If we take z = Zj in ( B.12 ) to be one of the turning points, then this equation al- 



lows us to determine the constants 4'o{zj). However, this gives only TV equations 
for the + 1 unknowns (jjQ^Zj), and also leaves the constant 7^0 undetermined. 
On a heuristic level, this means that in constructing the mapping one is free 
to choose both the origin and the scale. However, the sign of 7^0 is important, 
since it determines whether we are dealing with a barrier-type, or a well-type 
problem. Equation (B.12) also determines the mapping (t>Q{z) implicitly. 
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We can also choose Qi{(f>, h) to be a polynomial function, 



k=0 



Then equation (B.9) becomes 

M^) = l<t>'oRo'^' J^z' Ro'^' (^Ri ~ i<t>o?pKlk?j . (B.14) 

Now we require that this expression is non-singular at the turning points. This 
implies that [SH |32] 

dz' R-^'^ [r, ~ ^ 0g7fei ) = 0, (B.15) 



1/2 

y -111 - vv'o.' 

fe=0 / 

^(^oTfei I =0atz = zj, (B.16) 

fc=0 / 

where j runs from to iV and q runs from and rrij — 2. These conditions 
determine the constants 7^1. One should note that there are only /z— 1 equations, 
while there are constants. Therefore, some of them can be set to zero, and 
this considerably simplifies the expression for Qi{(j)). When these constants are 




determined, equation (B.14) determines (j)i{z). 

In a similar way, one can show that with j > 2 exist. Therefore, the 



mapping (B.7) is well-defined. Here we do not consider higher order corrections. 



for the purpose of the current paper it is enough to know that they exist [311 132] . 

Appendix B.2. Application to a first- order turning point 

In this subsection we illustrate the method explained in the previous section 
by the case of a single turning point on the real axis. We will use the method 
of comparison equations to solve the connection problem. We consider the case 
where the turning point separates a classically forbidden region (on the left) 



and a classically allowed region (on the right). From equation (10 1 for massless 
Dirac fermions, we see that 

Rq{x)^v'^{x)-pI, Ri{x) = iv'{x). (B.17) 

Since we consider a first-order turning point xq on the real axis, we have mp = 1 
and /i = 1, and for convenience we can set 0o(a;o) = 0. Since we assume that 
the classically forbidden region is on the left, we also have 710 = 1. Defining 
the square root as 

x^'^ = ^/^, X > 0, (B.18) 

a;i/2 = -iv^, a; < 0, (B.19) 
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from equation (B.12I we obtain 



(B.20) 



Since we have only one turning point, condition (B.15) does not give us any 
information, and we can set 711 — 701 = 0. Inserting these constants into 
equation (B.14|, we obtain (piix) as 



" Ri{x') 
Rl'\x') 



dx'. 



(B.21) 



Equations (B.20 1 and (B.21) together determine the mapping (l){x) 



Having estabhshed the mapping, we can now solve our related equation, 



-(t>V{(t)) = 0. 



(B.22) 



This is the well-known Airy (or Stokes) equation, and its solutions are given in 
terms of the Airy functions [4TJ |40j ; 

(B.23) 



V{(j)) = ciAi(-/i-2/30) + c2Bi(-/i-2/3^ 



Assuming that we are sufficiently far from the turning point, we can use the 
asymptotic expansions of the Airy functions. For ^ — >■ 00, they read 



Ai(e) 



_ 2 
3? 



e 3 



m) = 



63 - 



Ai(-0 = 



sin(§e 



2^3/2 , 



Bi(-0 = 



Cos(|e^/2+i^) 



(B.24) 
(B.25) 



We can then find the solutions to the original equation ( B.l ) from equation ( B.4 1 
and the mapping. From equation (B.7), we find that 



and from equation (B.20 1 that 



1/2 



(B.26) 



(B.27) 



We then construct the solution to the original equation using (B.4 1 and insert 
the mapping (B.20) and (B.21 1. Let us define 



1 



Rn 



Ri 



(B.28) 



Note that these waves are equal to the ones defined in equation (21). Further- 
more, they coincide with the basis functions (A.2) defined in Appendix A when 
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we do not regard the first exponent as an amplitude, but instead as part of 
the action, see the discussion in [Appendix A. 2 After a short computation, one 
obtains that for a: — ^ oo, 



t/4 



in I A: 



27r 



On the other hand, we find that for x 



27ri/2 



ciri^ix) 



-iir/4 



-1/2 



-C2?7i (a;), 



(B.29) 



(B.30) 



where Ro{x) in equation (B.28) is now negative and we have used the defini- 
tion (B.19I of the square root. 



In Appendix A| we stated that inside the classically forbidden region, only the 
term that increases along a given path should be kept [521 , and we put the other 
coefficient to zero by hand. Therefore one has to be careful when interpreting 



the results (B.29) and (B.30 1. When going from the classically allowed into the 



classically forbidden region, these equalities give [IHl [HI [23j 



where Cr 



= e"/'^7r-i/2(-ici + C2)/2 and c;,o 



Cl,oo)Vl 

= p-"/47r-l/2 



(B.31) 



TT-y'{lCi + C2)/2. 



Comparing this result with the ones obtained in Appendix A we conclude that 
the factor —i in front of Cr,oo on the right is nothing but the Stokes constant. 

Let us also see what happens when we go from the classically forbidden 
region into the classically allowed region. Then the connection formula reads 



?7 + (a;) r] + {x)+ir]^ (x). 



(B.32) 



This concludes our discussion of the connection formulae for a single turning 
point and their one-directional nature. 



Appendix B.3. Application to n-p and p-n junctions 

When we want to construct a uniform approximation for the transmission 
coefficient through an n-p junction, we should take into account two turning 
points X- < x+. We assume that they are nondegenerate, i.e. mg = mi — 1 
and fj, — 2. Since we consider a barrier type problem, we have 720 — 1, and we 
set (j)Q{x-) — — 0o(a^+) — ~o-- Defining the square root as in equation (A. 24), 
one finds from equation ( |B.12 ) that 



TTa 



\Jpy — v'^{x) dx — K, 



(B.33) 



where the last equality is implied by equation (61 ). This determines the constant 



a in the related equation in terms of the parameters of our initial problem. We 



use the same equation (B.12 1 to determine the mapping (t)o{z). Then we formally 
use the expansion 0o ^ ol, though the mapping may not be defined for such 
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^0- However, this does not influence the coeflicients in front of the asymptotic 
solutions. One finds that the mapping (f'oi^) is given by 



/v^x')-pldx'9,--cj,l + - + -\ni-2 



X-\- 



- ^ - ^ In 2 



X < X- 



X > x+. 



(B.34) 



For the case of an ordinary Schrodinger equation, one has i?i = and hence the 
first correction ipi is zero. This leads to the scattering matrices stated in [23] . 
However, for massless Dirac fermions (f>i does not vanish, since -Ri = iv'{x). 
Inserting this into equation (B.15I, and setting the upper limit to a;+, we find 
that 



701 = «• 



From equation (B.14), we then obtain (f>i{x) as 



01 (x) = 
01 (x) = 



1 

200 

1 

200 



: dx' + i In 



;2 - pi 



a 
200 



dx' — I In 



200 
a 



X < x^ 



X > x+, 



where we have once more made an expansion for large 0o. 
We find that the related equation (B.3l reduces to 



(B.35) 



(B.36) 



(B.37) 



which is exactly the equation for an n-p junction in graphene with a linear 
potential [11] [Ml [48] . Its solution is given by 

■ 2 

ViO = ciD,(V2e'"/4/i-i/20) + C2i?-,-i(v^e3^"/4/i-i/20), = ^, (B.38) 

where D^{x) are the parabolic cylinder functions [13 HI]. Its asymptotic ex- 
pansions are given by 



/27r 



2tt 



- 7r/2 < arg(z) < 7r/2 
arg (z) < — 7r/2 

arg (z) > 7r/2 



r(-^)' 

(B.39) 

Just as in the previous subsection, the solution of the original differential equa- 
tion is now given by 

^(x) = (0'(x))-i/V(0(x)), (B.40) 
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in which one has to insert the mapping (B.34), (B.36). Then we match the 



resulting expressions to the scattering states ( 30 ) . After some calculations, one 



finds that the transfer matrix ( 43 ) connecting the hole states on the right and 



the electron states on the left is given by 



T — 



le- 



i-n/2 



(B.41) 



where K is given by equation (61 1 and 9 by equation (73 1. To find the transfer 



matrix for a p-n junction, one can either do a similar calculation, or use the 



connection between the transfer matrices (43) and (47). Either way, one obtains 



T 



(B.42) 



From these matrices, one can easily derive the transmission coefficient (62) and 
the reflection coefficient (69 1. Finally, note that in the limit K/h — > oo, the 



equation (A. 28) 



matrix (B.41) reduces to the matrices obtained from the WKB approximation. 



Appendix B.4- Application to Schrddinger-like cases 

In the previous subsection, one of the turning points corresponded to v{z) = 
— \py\, and the other one to v{z) = \py\. In this subsection we consider the 
situation where the two turning points correspond to v{z) — — \py \ . This includes 
both the conventional tunneling regime from section [6j where we took only the 
two real turning points into account, as well as the regime of above-barrier 
scattering from section [5] when we consider only the middle two turning points. 
The final answers for the refiection and transmission coefficients in this case are 
similar to those for an ordinary Schrodinger equation. Therefore we speak of 
Schrodinger-like cases. 

The computations for both cases are similar P3|, but the one for above- 
barrier scattering is slightly more complicated. Therefore we will focus on it 
below, leaving the other one to the reader. For above-barrier scattering the 
turning points are complex. Hence we have to consider the method of compar- 
ison equations in a complex domain D, containing these turning points as well 
as the real axis. 

Let us consider two simple complex turning points, zi_ and Since the 
problem is of barrier-type, 720 = 1, and we set 0o(-2i-) = ^'P{zi+) = We 
assume that the branch cut is placed between the two turning points, and define 

-r.1/2 = ^fz 



rl/2 



/X, 



X > Xq, 
X < Xq, 



(B.43) 
(B.44) 



where xq is the point where the Stokes line from to zi_ crosses the real 
axis, see section [5] Performing the integration on the positive side of the cut, 
we find that 



K ^ 2i 



21+ 



{v\z)~plf'^dz = i {v\z) 



21 + 



-Plf'dz^ 



(B.45) 
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We remind the reader that xq is thepoint where the Stokes hne from zi_|_ to 
zi- crosses the real axis, see section pj and that with x^, we mean the point 
Xo zt e when e — > 0. The first equahty above fohows from the definition in 



equation (79), and the second one from the fact that v{x) is a real function. 



The third equality, which is the most important for this subsection, follows 



from equation (B.12I 



With the help of equation (B.12|, one then determines the mapping (j^oiz)- 
We set the lower limit to zi-, and then split the integrals into two parts. The 
first part goes along the cut, and connects the turning point to the real axis. 
Using equation (B.45) and expanding for large 0o, we find that 



v^{xr)-pldx'^-l^l-'--';ini-2^' 



I 1 



4 2 

,2 52 



2 \ b 



X < Xo 



X > Xo- 



(B.46) 



Contrary to the two previous examples, equation (B.15I shows that this time 



701 = 0, 



(B.47) 



due to the fact that both turning points correspond to v{zi±) — —\py\. This 
implies that the comparison equation will be identical with the one for an or- 
dinary Schrodinger equation, and that we can therefore expect similar results. 
Using equation (B.14), we find that 0i(a;) equals 



h{x) = -\^o^ 



v{x) 
-\Pv\ 



idv 



— Py 



v{x) 



\Pv\ 



idv 



X < Xo, 



X > Xo- 



(B.48) 



The related equation (B.3) reduces to 



,d^V 



(B.49) 



which is indeed the same as for an ordinary Schrodinger equation. Its solution 
is given by 



ib^ 



(B.50) 

As in the previous subsection, we now construct the exact solution using equa- 
tion (B.4). Then we make asymptotic expansions of the parabolic cylinder 
functions and apply the mapping (B.46 1, (B.48). Matching the result to the 



scattering states ( 30 1, we obtain the reflection and transmission coefficients ( 82 ) 



which are similar to those for an ordinary Schrodinger equation |23j . 
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The computation for tunneling through the barrier runs entirely similar, the 
main difference being that this time the turning points are real and that K 
is positive instead of negative. However, as indicated in section [6] and shown 
in [33], up to these differences the final answer is exactly the same. 



Appendix B.5. Application to above-barrier scattering 

In the previous subsection we applied the method of comparison equations 
to above-barrier scattering for two turning points that both correspond to 
v{z) = — \py\- However, in section[5j we saw that upon near-normal incidence on 
a short-range potential, the two complex turning points zi+ and Z2+, that cor- 
respond to v{zi+) = —\py\ and v{z2+) = \py\ respectively, merge. So if we want 
to derive an expression for the reflection coefficient that is valid at near-normal 
incidence, we should apply the method of comparison equations to these two 
turning points. Since they are connected by an anti-Stokes line, we are dealing 
with a well- type problem, and hence 720 = — 1. The turning points and Z2+ 
are mapped to —ib and ib respectively, and the branch cut is placed between 



these two points. Applying the mapping (B.12), we find that 



5 = 



iv\z)-piy/^dz = -'^.. 



(B.51) 



where S was defined earlier in equation (87). From figure [3|b), we see that 
four anti-Stokes lines emerge from the cluster. Since we are interested in the 
wave function along the lower two lines, we consider 0o in the lower half-plane. 



Applying equation (B.12 1 once more, we find 



^1+ 



iv' - piy/'dz' = ~ + \M^l + b^f' + ^ In 



0| 
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We proceed by determining 701 with equation (B.15I, and obtain 701 — 1. This 
allows us to calculate (piiz). We find that 



1/2- 



(B.52) 




iv'{z') 



dz' + i In 







\ 1/2' 


-5) 


00 






b 


^62 





The related equation (|B.3) reduces to 

+ (-02 -b'^ + h) V{0 = 0, 



,d''V 



(B.53) 



(B.54) 



and has the solution 

V{(l)) = CiD^{V2h'^/^q^) + C2C-.-i(zV2/i-i/2 



2h' 



(B.55) 



in terms of the aforementioned parabolic cylinder functions. Looking at the 



asymptotic expansions (B.39), we see the solutions represent traveling waves 
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along the lines Arg((/)) = — 7r/4 and Arg((^) — — 37r/4 in the lower half plane. 
Therefore we make an expansion along these lines, and then make use of the 



mapping (B.52) and (B.53I. We then obtain the asymptotic expansion of the 



solution of the original equation on the anti -Sto kes lines 71 and 72 in figure [10} 
In terms of the functions rj^{z), defined in (21 1, let us write 



on 7i and 



on 72- Then we find that 



(B.56) 



(B.57) 



(B.58) 



where a was defined in equation (86). The first exponent in (21 1 becomes 
{—g)^^^"^ , and we immediately see that upon passing from 72 to 71 , equation ( 84 1 



goes over in equation ( 85 ) 



Appendix B.6. Application to the exactly solvable potential increase 

In this final subsection, we apply the method of comparison equations to the 
situation considered in section[7] As before, we consider a cluster of two turning 
points, namely the turning points Z2+ and Z14. that lie above the real axis, 
with 'Re[z2+) < Re(zn_), and v{z2+) = \py\, v{zi^) — ~\py\- These turning 
points are mapped to —a and a respectively, and we see that the situation 
exactly coincides with the p-n junction that was considered in [Appendix B.3[ 
However, this time we want to consider the asymptotic representation of the 
exact solutions along different anti-Stokes lines. To find the transmitted wave, 
we need the solution along the anti-Stokes line 71 in figure[TT] which corresponds 
to Arg((^) = in the comparison equation. On the other hand, we find the 
incoming and reflected waves by considering the solution along the anti-Stokes 
line 72, that corresponds to the line Arg(0) = —tt/2 in the comparison equation. 
A second difference is that we want to have all wavefunctions defined with 
respect to the reference point zi+. 

Applying the method of comparison equations as in Appendix B.3 one can 
show that the comparison equation for this case equals 



which has the solution 



2h 



(B.59) 



, (B.60) 



in terms of the parabolic cylinder functions. We then use the asymptotic ex- 
pansions of the parabolic cylinder functions and apply the mapping — (j){z) to 
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find the asymptotic expansion of the exact solution along the anti-Stokes lines 
Introducing the coefficients cj^ and cl^2 the previous 
'that 



7i and 72 in figure 



11 



subsection, one finos 



„7l 



„7l 



-2S/h 



1/2 



-i0 



(B.61) 



in terms of the quantities S and 6 defined by equations (117) and ( 118 ), respect- 
ively. 



The reflection coefficient (116) is then obtained by a procedure similar to 



that in section [5j One starts with a transmitted wave, that is defined with 
respect to the point xa on the real axis. Then one changes the reference point 



to and uses the result (B.61) to make the transition from the anti-Stokes 



line 7i to the anti-Stokes line 72. Finally, one changes the reference point back 



to xq and finds the refiection coefficient (116) 



Appendix C. Relations betv^reen asymptotic scattering states 



In a previous paper 14 the authors placed particular emphasis on the geo- 
metric interpretation of the amplitude factor of the function ^'(a;) in (|4| in terms 
of the Berry phase. The scattering states defined there are different from those 
we have introduced in the current paper. Here we establish the relationship 
between different scattering states. 

First we notice that using (11 1 we can write relation ([9| as 



* = 



«sgn(py).g 



1 

\Py\ 



dx 



v(x) 



Using equation ( 30 ) and ( 34 1 we find 



[G^'/\x)~t^sgn{py)G^^/Hx) 
In the electron region (i/ = — 1) the latter can be written as 



where (j)'^ 



^± = 

j)- = 7rsgn(p.y), 
cos(0+) 



\Py\ 



Px 

v{x)y 



Py 

Hx)[ 



For the hole region {v = 1) we can write (C.2) as 



p±iS(^xo,x)/h 
\]>_l_ = p'^sgn(p„)/2 



\/Px{x) 



G^^/^{x)-isgn{py)G^^'^{x) 
-G^^/^{x) - isgn{py)G^^/^{x) 



(C.l) 



(C.2) 



(C.3) 



(C.4) 



(C.5) 
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or 



„±iS{x,xo)/h \2\v(t]\ ( p-'"^J/2 \ 



Comparing (C.3), (C.6) with asymptotic scattering states defined in [T3] we 



find that in the electron region 

^± = ^\Vy\-^'^^±, (C.7) 

while in the hole region 

The change in the definition of the asymptotic scattering states leads to a cor- 
responding change of the phase of the transmission coefficient for tunneling from 



an electron region to a hole region, cf . equation ( 62 ) of the current paper and 
equation (125) in [14^ . 
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